(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 9.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[    165247,       4341]
NotebookOptionsPosition[    158610,       4142]
NotebookOutlinePosition[    158975,       4158]
CellTagsIndexPosition[    158932,       4155]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Setup MEH[k]\[LineSeparator]", "Subsubsection",
 CellChangeTimes->{{3.570327517917389*^9, 3.57032751997117*^9}, {
   3.570327554725595*^9, 3.57032755584058*^9}, {3.570327787927311*^9, 
   3.5703277888672857`*^9}, {3.5703722949756393`*^9, 3.570372296323347*^9}, {
   3.592769661830781*^9, 3.592769666954217*^9}, {3.5934466262643137`*^9, 
   3.593446628768866*^9}, {3.593469393516514*^9, 3.593469394584261*^9}, 
   3.593525800106504*^9, {3.593964656104784*^9, 3.593964657792479*^9}, {
   3.5944913560479097`*^9, 3.594491356878201*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"c", " ", "=", " ", "299792.458"}], ";", " ", 
  RowBox[{"z", " ", "=", " ", "0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"h", "=", "0.70"}], ";", 
  RowBox[{"Omm", " ", "=", " ", "0.27"}], ";", 
  RowBox[{"Omb", "=", "0.04"}], ";", 
  RowBox[{"keq", " ", "=", " ", 
   FractionBox[
    RowBox[{"0.0746", "*", "Omm", "*", 
     SuperscriptBox["h", "2"]}], 
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{"2.728", "/", "2.7"}], ")"}], "2"]]}], ";", 
  RowBox[{"zdrag", "=", 
   RowBox[{"1291", "*", 
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{"Omm", "*", 
        RowBox[{"h", "^", "2"}]}], ")"}], "^", 
      RowBox[{"(", "0.251", ")"}]}], "/", 
     RowBox[{"(", 
      RowBox[{"1", "+", 
       RowBox[{"0.659", "*", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"Omm", "*", 
           RowBox[{"h", "^", "2"}]}], ")"}], "^", 
         RowBox[{"(", "0.828", ")"}]}]}]}], ")"}]}], "*", 
    RowBox[{"(", 
     RowBox[{"1", "+", 
      RowBox[{"0.313", "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"Omm", "*", 
          RowBox[{"h", "^", "2"}]}], ")"}], "^", 
        RowBox[{"(", 
         RowBox[{"-", "0.419"}], ")"}]}], "*", 
       RowBox[{"(", 
        RowBox[{"1", "+", 
         RowBox[{"0.607", "*", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"Omm", "*", 
             RowBox[{"h", "^", "2"}]}], ")"}], "^", "0.674"}]}]}], ")"}], "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"Omb", "*", 
          RowBox[{"h", "^", "2"}]}], ")"}], "^", 
        RowBox[{"(", 
         RowBox[{"0.238", "*", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"Omm", "*", 
             RowBox[{"h", "^", "2"}]}], ")"}], "^", "0.223"}]}], ")"}]}]}]}], 
     ")"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"zeq", "=", 
   RowBox[{"2.50", "*", 
    RowBox[{"10", "^", "4"}], "*", "Omm", "*", 
    RowBox[{
     RowBox[{"h", "^", "2"}], "/", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"2.726", "/", "2.7"}], ")"}], "^", "4"}]}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"Rd", "=", 
    RowBox[{"31.5", "*", "Omb", "*", 
     RowBox[{"h", "^", "2"}], "*", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"2.726", "/", "2.7"}], ")"}], "^", 
       RowBox[{"(", 
        RowBox[{"-", "4"}], ")"}]}], "/", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"1", "+", "zdrag"}], ")"}], "/", 
        RowBox[{"10", "^", "3"}]}], ")"}]}]}]}], ";"}], 
  RowBox[{"(*", " ", 
   RowBox[{"R_drag", " ", "in", " ", 
    RowBox[{"tf_fit", ".", "c"}]}], " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Req", "=", 
   RowBox[{"31.5", "*", "Omb", "*", 
    RowBox[{"h", "^", "2"}], "*", 
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{"2.726", "/", "2.7"}], ")"}], "^", 
      RowBox[{"(", 
       RowBox[{"-", "4"}], ")"}]}], "/", 
     RowBox[{"(", 
      RowBox[{"zeq", "/", 
       RowBox[{"10", "^", "3"}]}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", 

 RowBox[{
  RowBox[{
   RowBox[{"kSilk", "=", 
    RowBox[{"1.6", "*", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"Omb", "*", 
        RowBox[{"h", "^", "2"}]}], ")"}], "^", 
      RowBox[{"(", "0.52", ")"}]}], "*", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"Omm", "*", 
        RowBox[{"h", "^", "2"}]}], ")"}], "^", 
      RowBox[{"(", "0.73", ")"}]}], "*", 
     RowBox[{"(", 
      RowBox[{"1", "+", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"10.4", "*", "Omm", "*", 
          RowBox[{"h", "^", "2"}]}], ")"}], "^", 
        RowBox[{"(", 
         RowBox[{"-", "0.95"}], ")"}]}]}], ")"}]}]}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{"in", " ", 
    RowBox[{"Mpc", "^", 
     RowBox[{"(", 
      RowBox[{"-", "1"}], ")"}]}]}], " ", "*)"}]}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"xx", "=", 
    RowBox[{
     RowBox[{
      RowBox[{"2", "/", "3"}], "/", "keq"}], "*", 
     RowBox[{"Sqrt", "[", 
      RowBox[{"6", "/", "Req"}], "]"}], "*", 
     RowBox[{"Log", "[", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"Sqrt", "[", 
          RowBox[{"1", "+", "Rd"}], "]"}], "+", 
         RowBox[{"Sqrt", "[", 
          RowBox[{"Rd", "+", "Req"}], "]"}]}], ")"}], "/", 
       RowBox[{"(", 
        RowBox[{"1", "+", 
         RowBox[{"Sqrt", "[", "Req", "]"}]}], ")"}]}], "]"}]}]}], ";"}], 
  StyleBox[
   RowBox[{"(*", " ", "sound_horizon", " ", "*)"}],
   FontColor->GrayLevel[0]]}], "\n", 
 RowBox[{
  RowBox[{"\[Beta]c", "=", 
   RowBox[{"1", "/", 
    RowBox[{"(", 
     RowBox[{"1", "+", 
      RowBox[{
       RowBox[{"0.944", "/", 
        RowBox[{"(", 
         RowBox[{"1", "+", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"458", "*", "Omm", "*", 
             RowBox[{"h", "^", "2"}]}], ")"}], "^", 
           RowBox[{"(", 
            RowBox[{"-", "0.708"}], ")"}]}]}], ")"}]}], "*", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"Omm", "-", "Omb"}], ")"}], "/", "Omm"}], ")"}], "^", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"0.395", "*", "Omm", "*", 
              RowBox[{"h", "^", "2"}]}], ")"}], "^", 
            RowBox[{"(", 
             RowBox[{"-", "0.0266"}], ")"}]}], ")"}]}], "-", "1"}], ")"}]}]}],
      ")"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Beta]b", "=", 
   RowBox[{"0.5", "+", 
    RowBox[{"Omb", "/", "Omm"}], "+", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"3", "-", 
       RowBox[{"2", "*", 
        RowBox[{"Omb", "/", "Omm"}]}]}], ")"}], "*", 
     RowBox[{"Sqrt", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"17.2", "*", "Omm", "*", 
          RowBox[{"h", "^", "2"}]}], ")"}], "^", "2"}], "+", "1"}], 
      "]"}]}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Alpha]c", "=", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"46.9", "*", "Omm", "*", 
          RowBox[{"h", "^", "2"}]}], ")"}], "^", 
        RowBox[{"(", "0.670", ")"}]}], "*", 
       RowBox[{"(", 
        RowBox[{"1", "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"32.1", "*", "Omm", "*", 
            RowBox[{"h", "^", "2"}]}], ")"}], "^", 
          RowBox[{"(", 
           RowBox[{"-", "0.532"}], ")"}]}]}], ")"}]}], ")"}], "^", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "Omb"}], "/", "Omm"}], ")"}]}], "*", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"12.0", "*", "Omm", "*", 
          RowBox[{"h", "^", "2"}]}], ")"}], "^", 
        RowBox[{"(", "0.424", ")"}]}], "*", 
       RowBox[{"(", 
        RowBox[{"1", "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"45.0", "*", "Omm", "*", 
            RowBox[{"h", "^", "2"}]}], ")"}], "^", 
          RowBox[{"(", 
           RowBox[{"-", "0.582"}], ")"}]}]}], ")"}]}], ")"}], "^", 
     RowBox[{"(", 
      RowBox[{"-", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"Omb", "/", "Omm"}], ")"}], "^", "3"}]}], ")"}]}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Alpha]b", "=", 
   RowBox[{"2.07", "*", "keq", "*", "xx", "*", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"1", "+", "Rd"}], ")"}], "^", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "3"}], "/", "4"}], ")"}]}], "*", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"1", "+", "zeq"}], ")"}], "/", 
     RowBox[{"(", 
      RowBox[{"1", "+", "zdrag"}], ")"}]}], "*", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       RowBox[{"-", "6"}], "*", 
       RowBox[{"Sqrt", "[", 
        RowBox[{"1", "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"1", "+", "zeq"}], ")"}], "/", 
          RowBox[{"(", 
           RowBox[{"1", "+", "zdrag"}], ")"}]}]}], "]"}]}], "+", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"2", "+", 
         RowBox[{"3", "*", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"1", "+", "zeq"}], ")"}], "/", 
           RowBox[{"(", 
            RowBox[{"1", "+", "zdrag"}], ")"}]}]}]}], ")"}], "*", 
       RowBox[{"Log", "[", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"Sqrt", "[", 
            RowBox[{"1", "+", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{"1", "+", "zeq"}], ")"}], "/", 
              RowBox[{"(", 
               RowBox[{"1", "+", "zdrag"}], ")"}]}]}], "]"}], "+", "1"}], 
          ")"}], "/", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"Sqrt", "[", 
            RowBox[{"1", "+", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{"1", "+", "zeq"}], ")"}], "/", 
              RowBox[{"(", 
               RowBox[{"1", "+", "zdrag"}], ")"}]}]}], "]"}], "-", "1"}], 
          ")"}]}], "]"}]}]}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"TransferEH", "[", "k_", "]"}], ":=", " ", 
   RowBox[{"Piecewise", "[", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{"Omm", "-", "Omb"}], ")"}], "/", "Omm"}], "*", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{"1", "/", 
               RowBox[{"(", 
                RowBox[{"1", "+", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{"k", "*", 
                    RowBox[{"xx", "/", "5.4"}]}], ")"}], "^", "4"}]}], 
                ")"}]}], ")"}], "*", 
             RowBox[{
              RowBox[{"Log", "[", 
               RowBox[{"E", "+", 
                RowBox[{"1.8", "*", "\[Beta]c", "*", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"k", "/", "keq"}], "/", "13.41"}], ")"}]}]}], 
               "]"}], "/", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{"Log", "[", 
                 RowBox[{"E", "+", 
                  RowBox[{"1.8", "*", "\[Beta]c", "*", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"k", "/", "keq"}], "/", "13.41"}], ")"}]}]}], 
                 "]"}], "+", 
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{"14.2", "+", 
                   RowBox[{"386", "/", 
                    RowBox[{"(", 
                    RowBox[{"1", "+", 
                    RowBox[{"69.9", "*", 
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"k", "/", "keq"}], "/", "13.41"}], ")"}], "^", 
                    "1.08"}]}]}], ")"}]}]}], ")"}], "*", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"k", "/", "keq"}], "/", "13.41"}], ")"}], "^", 
                  "2"}]}]}], ")"}]}]}], "+", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"1", "-", 
               RowBox[{"(", 
                RowBox[{"1", "/", 
                 RowBox[{"(", 
                  RowBox[{"1", "+", 
                   RowBox[{
                    RowBox[{"(", 
                    RowBox[{"k", "*", 
                    RowBox[{"xx", "/", "5.4"}]}], ")"}], "^", "4"}]}], 
                  ")"}]}], ")"}]}], ")"}], "*", 
             RowBox[{
              RowBox[{"Log", "[", 
               RowBox[{"E", "+", 
                RowBox[{"1.8", "*", "\[Beta]c", "*", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"k", "/", "keq"}], "/", "13.41"}], ")"}]}]}], 
               "]"}], "/", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{"Log", "[", 
                 RowBox[{"E", "+", 
                  RowBox[{"1.8", "*", "\[Beta]c", "*", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"k", "/", "keq"}], "/", "13.41"}], ")"}]}]}], 
                 "]"}], "+", 
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"14.2", "/", "\[Alpha]c"}], "+", 
                   RowBox[{"386", "/", 
                    RowBox[{"(", 
                    RowBox[{"1", "+", 
                    RowBox[{"69.9", "*", 
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"k", "/", "keq"}], "/", "13.41"}], ")"}], "^", 
                    "1.08"}]}]}], ")"}]}]}], ")"}], "*", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"k", "/", "keq"}], "/", "13.41"}], ")"}], "^", 
                  "2"}]}]}], ")"}]}]}]}], ")"}]}], "+", 
         RowBox[{
          RowBox[{"Omb", "/", "Omm"}], "*", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"Log", "[", 
                RowBox[{"E", "+", 
                 RowBox[{"1.8", "*", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"k", "/", "keq"}], "/", "13.41"}], ")"}]}]}], 
                "]"}], "/", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{"Log", "[", 
                  RowBox[{"E", "+", 
                   RowBox[{"1.8", "*", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"k", "/", "keq"}], "/", "13.41"}], ")"}]}]}], 
                  "]"}], "+", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{"14.2", "+", 
                    RowBox[{"386", "/", 
                    RowBox[{"(", 
                    RowBox[{"1", "+", 
                    RowBox[{"69.9", "*", 
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"k", "/", "keq"}], "/", "13.41"}], ")"}], "^", 
                    RowBox[{"(", "1.08", ")"}]}]}]}], ")"}]}]}], ")"}], "*", 
                  RowBox[{
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"k", "/", "keq"}], "/", "13.41"}], ")"}], "^", 
                   "2"}]}]}], ")"}]}], ")"}], "/", 
             RowBox[{"(", 
              RowBox[{"1", "+", 
               RowBox[{
                RowBox[{"(", 
                 RowBox[{"k", "*", 
                  RowBox[{"xx", "/", "5.2"}]}], ")"}], "^", "2"}]}], ")"}]}], 
            "+", 
            RowBox[{"\[Alpha]b", "*", 
             RowBox[{
              RowBox[{"Exp", "[", 
               RowBox[{"-", 
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{"k", "/", "kSilk"}], ")"}], "^", "1.4"}]}], "]"}], 
              "/", 
              RowBox[{"(", 
               RowBox[{"1", "+", 
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"\[Beta]b", "/", "k"}], "/", "xx"}], ")"}], "^", 
                 "3"}]}], ")"}]}]}]}], ")"}], 
          RowBox[{
           RowBox[{"Sin", "[", 
            RowBox[{"k", "*", 
             RowBox[{"xx", "/", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{"1", "+", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{"8.41", "*", 
                    RowBox[{
                    RowBox[{
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{"Omm", "*", 
                    RowBox[{"h", "^", "2"}]}], ")"}], "^", 
                    RowBox[{"(", "0.435", ")"}]}], "/", "k"}], "/", "xx"}]}], 
                   ")"}], "^", "3"}]}], ")"}], "^", 
               RowBox[{"(", 
                RowBox[{"1", "/", "3"}], ")"}]}]}]}], "]"}], "/", 
           RowBox[{"(", 
            RowBox[{"k", "*", 
             RowBox[{"xx", "/", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{"1", "+", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{"8.41", "*", 
                    RowBox[{
                    RowBox[{
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{"Omm", "*", 
                    RowBox[{"h", "^", "2"}]}], ")"}], "^", 
                    RowBox[{"(", "0.435", ")"}]}], "/", "k"}], "/", "xx"}]}], 
                   ")"}], "^", "3"}]}], ")"}], "^", 
               RowBox[{"(", 
                RowBox[{"1", "/", "3"}], ")"}]}]}]}], ")"}]}]}]}], ",", 
        RowBox[{"k", ">", "0"}]}], "}"}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"1.0", ",", " ", 
        RowBox[{"k", "\[Equal]", "0"}]}], "}"}]}], "}"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"Hubble", "[", "zz_", "]"}], ":=", 
    RowBox[{"100", " ", "h", " ", 
     SqrtBox[
      RowBox[{
       RowBox[{"Omm", " ", 
        SuperscriptBox[
         RowBox[{"(", 
          RowBox[{"1", "+", "zz"}], ")"}], "3"]}], "+", 
       RowBox[{"(", 
        RowBox[{"1", "-", "Omm"}], ")"}]}]]}]}], ";"}], "\n", 
  "\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"don", "'"}], "t", " ", "forget", " ", "below", " ", "is", " ", 
    "only", " ", "valid", " ", "for", " ", "a", " ", "flat", " ", 
    "\[CapitalLambda]CDM", " ", "or", " ", "CDM", " ", "universe"}], 
   "*)"}]}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"Growth", "[", "zz_", "]"}], ":=", 
   RowBox[{
    FractionBox[
     RowBox[{"5", " ", "Omm"}], "2"], " ", 
    FractionBox[
     RowBox[{"Hubble", "[", "zz", "]"}], 
     RowBox[{"100", " ", "h"}]], 
    RowBox[{"NIntegrate", "[", 
     RowBox[{
      FractionBox[
       RowBox[{"(", 
        RowBox[{"1", "+", "y"}], ")"}], 
       SuperscriptBox[
        RowBox[{"(", 
         FractionBox[
          RowBox[{"Hubble", "[", "y", "]"}], 
          RowBox[{"100", " ", "h"}]], ")"}], "3"]], ",", 
      RowBox[{"{", 
       RowBox[{"y", ",", "zz", ",", "\[Infinity]"}], "}"}]}], "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"MEH", "[", "k_", "]"}], ":=", 
    RowBox[{
     RowBox[{
      RowBox[{"2.", "/", "3."}], "/", "Omm"}], "*", 
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{"c", "/", 
        RowBox[{"(", 
         RowBox[{"100", " ", "h"}], ")"}]}], ")"}], "2"], "*", 
     RowBox[{"TransferEH", "[", "k", "]"}], "*", 
     SuperscriptBox["k", "2"], "*", 
     RowBox[{"Growth", "[", "z", "]"}]}]}], ";"}], "\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{"Simple", " ", 
    RowBox[{"Transfer", "[", "k", "]"}]}], " ", 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Omb", " ", "=", " ", "0.04"}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{"baryon", " ", "density"}], " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[CapitalGamma]", "=", 
   RowBox[{"Omm", "*", "h", "*", 
    RowBox[{"Exp", "[", 
     RowBox[{
      RowBox[{"-", "Omb"}], 
      RowBox[{"(", 
       RowBox[{"1", "+", 
        RowBox[{
         RowBox[{"Sqrt", "[", 
          RowBox[{"2", "*", "h"}], "]"}], "/", "Omm"}]}], ")"}]}], "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Transfer", "[", "k_", "]"}], ":=", 
   RowBox[{
    RowBox[{
     RowBox[{"Log", "[", 
      RowBox[{"1", "+", 
       RowBox[{"2.34", 
        RowBox[{"(", 
         RowBox[{"k", "/", 
          RowBox[{"(", 
           RowBox[{"h", "*", "\[CapitalGamma]"}], ")"}]}], ")"}]}]}], "]"}], 
     "/", 
     RowBox[{"(", 
      RowBox[{"2.34", "*", 
       RowBox[{"(", 
        RowBox[{"k", "/", 
         RowBox[{"(", 
          RowBox[{"h", "*", "\[CapitalGamma]"}], ")"}]}], ")"}]}], ")"}]}], 
    "*", 
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{"1", "+", 
       RowBox[{"3.89", "*", 
        RowBox[{"(", 
         RowBox[{"k", "/", 
          RowBox[{"(", 
           RowBox[{"h", "*", "\[CapitalGamma]"}], ")"}]}], ")"}]}], "+", 
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{"16.1", "*", 
          RowBox[{"(", 
           RowBox[{"k", "/", 
            RowBox[{"(", 
             RowBox[{"h", "*", "\[CapitalGamma]"}], ")"}]}], ")"}]}], ")"}], 
        "2"], "+", 
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{"5.46", "*", 
          RowBox[{"(", 
           RowBox[{"k", "/", 
            RowBox[{"(", 
             RowBox[{"h", "*", "\[CapitalGamma]"}], ")"}]}], ")"}]}], ")"}], 
        "3"], "+", 
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{"6.71", "*", 
          RowBox[{"(", 
           RowBox[{"k", "/", 
            RowBox[{"(", 
             RowBox[{"h", "*", "\[CapitalGamma]"}], ")"}]}], ")"}]}], ")"}], 
        "4"]}], ")"}], 
     RowBox[{"-", "0.25"}]]}]}], ";"}]}], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.563542383318079*^9, 3.563542460983938*^9}, 
   3.563542493723468*^9, {3.563542586887805*^9, 3.563542628749124*^9}, {
   3.563543087042822*^9, 3.563543090748076*^9}, {3.563543374436173*^9, 
   3.563543422198125*^9}, {3.5644021704266863`*^9, 3.564402172967256*^9}, {
   3.564405767539486*^9, 3.564405768052952*^9}, {3.564413235949687*^9, 
   3.564413236853958*^9}, {3.564682427610278*^9, 3.564682436448017*^9}, {
   3.5646824745485983`*^9, 3.564682480322901*^9}, {3.5648343000988703`*^9, 
   3.56483430413787*^9}, {3.564850105889791*^9, 3.56485010797779*^9}, {
   3.566426830196505*^9, 3.566426831159525*^9}, 3.570327537517303*^9, {
   3.570327712775303*^9, 3.570327714117038*^9}, {3.5703768187944307`*^9, 
   3.570376821241983*^9}, {3.5715755872881193`*^9, 3.571575636954029*^9}, {
   3.5885027348150806`*^9, 3.588502735071302*^9}, {3.588502813821426*^9, 
   3.588502814028502*^9}, {3.5885032878929996`*^9, 3.588503339433796*^9}, {
   3.588503406456326*^9, 3.588503406639093*^9}, {3.588503465987542*^9, 
   3.588503471923272*^9}, {3.5922255614844513`*^9, 3.592225562129325*^9}, {
   3.5927643164655037`*^9, 3.592764317763259*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"MEH", "[", 
  RowBox[{"0.01", "*", "h"}], "]"}]], "Input",
 CellChangeTimes->{{3.592152851087934*^9, 3.592152852570692*^9}, {
  3.592152889472789*^9, 3.592152895286504*^9}}],

Cell[BoxData["1275.1632126042107`"], "Output",
 CellChangeTimes->{3.592152895936891*^9, 3.592223409910748*^9, 
  3.592772110419016*^9, 3.592772234865951*^9, 3.5929147545551453`*^9, 
  3.5931914766945868`*^9, 3.593385026671564*^9, 3.593387544136036*^9, 
  3.593443816729603*^9, 3.593443973908484*^9, 3.593445028950808*^9, 
  3.593446173308908*^9, 3.593446638065283*^9, 3.593469399136197*^9, 
  3.593525811661569*^9, 3.5935259825637693`*^9, 3.593963581329606*^9, 
  3.593964662793477*^9, 3.594040331041827*^9, 3.5940408505378017`*^9, 
  3.59404304436619*^9, 3.5944913624800043`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"LogLogPlot", "[", 
  RowBox[{
   RowBox[{
    SuperscriptBox[
     RowBox[{"MEH", "[", 
      RowBox[{"k", "*", "h"}], "]"}], "2"], "*", 
    RowBox[{"Pvar", "[", "k", "]"}]}], ",", " ", 
   RowBox[{"{", 
    RowBox[{"k", ",", " ", "0.001", ",", "1"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.592223642355777*^9, 3.592223770239456*^9}, {
   3.5922238479599323`*^9, 3.5922238482505207`*^9}, {3.592223900806583*^9, 
   3.592223902446669*^9}, 3.592223973952518*^9}],

Cell[BoxData[
 GraphicsBox[{},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesLabel->{None, None},
  AxesOrigin->{-6.907755278982137, -0.6931471805599453},
  CoordinatesToolOptions:>{"DisplayFunction" -> ({
      Exp[
       Part[#, 1]], 
      Exp[
       Part[#, 2]]}& ), "CopiedValueFunction" -> ({
      Exp[
       Part[#, 1]], 
      Exp[
       Part[#, 2]]}& )},
  FrameTicks->{{{{0., 
       FormBox[
        TagBox[
         InterpretationBox["\"1.0\"", 1., AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {2.302585092994046, 
       FormBox[
        TagBox[
         InterpretationBox["\"10.0\"", 10., AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {1.6094379124341003`, 
       FormBox[
        TagBox[
         InterpretationBox["\"5.0\"", 5., AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {3.912023005428146, 
       FormBox[
        TagBox[
         InterpretationBox["\"50.0\"", 50., AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {0.6931471805599453, 
       FormBox[
        TagBox[
         InterpretationBox["\"2.0\"", 2., AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {2.995732273553991, 
       FormBox[
        TagBox[
         InterpretationBox["\"20.0\"", 20., AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {1.0986122886681098`, 
       FormBox[
        TagBox[
         InterpretationBox["\"3.0\"", 3., AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {3.4011973816621555`, 
       FormBox[
        TagBox[
         InterpretationBox["\"30.0\"", 30., AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {0.4054651081081644, 
       FormBox[
        TagBox[
         InterpretationBox["\"1.5\"", 1.5, AutoDelete -> True], NumberForm[#, {
           DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {2.70805020110221, 
       FormBox[
        TagBox[
         InterpretationBox["\"15.0\"", 15., AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {1.9459101490553132`, 
       FormBox[
        TagBox[
         InterpretationBox["\"7.0\"", 7., AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {4.248495242049359, 
       FormBox[
        TagBox[
         InterpretationBox["\"70.0\"", 70., AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {0.09531017980432493, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.1823215567939548, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.26236426446749106`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.336472236621213, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.47000362924573563`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.5306282510621705, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.5877866649021191, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.6418538861723948, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {1.3862943611198906`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {1.791759469228055, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.0794415416798357`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.1972245773362196`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.3978952727983707`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.4849066497880004`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.5649493574615367`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.6390573296152584`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.772588722239781, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.833213344056216, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.8903717578961645`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.9444389791664403`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {3.6888794541139363`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {4.0943445622221, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}}, {{0., 
       FormBox["\"\"", TraditionalForm]}, {2.302585092994046, 
       FormBox["\"\"", TraditionalForm]}, {1.6094379124341003`, 
       FormBox["\"\"", TraditionalForm]}, {3.912023005428146, 
       FormBox["\"\"", TraditionalForm]}, {0.6931471805599453, 
       FormBox["\"\"", TraditionalForm]}, {2.995732273553991, 
       FormBox["\"\"", TraditionalForm]}, {1.0986122886681098`, 
       FormBox["\"\"", TraditionalForm]}, {3.4011973816621555`, 
       FormBox["\"\"", TraditionalForm]}, {0.4054651081081644, 
       FormBox["\"\"", TraditionalForm]}, {2.70805020110221, 
       FormBox["\"\"", TraditionalForm]}, {1.9459101490553132`, 
       FormBox["\"\"", TraditionalForm]}, {4.248495242049359, 
       FormBox["\"\"", TraditionalForm]}, {0.09531017980432493, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.1823215567939548, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.26236426446749106`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.336472236621213, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.47000362924573563`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.5306282510621705, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.5877866649021191, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.6418538861723948, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {1.3862943611198906`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {1.791759469228055, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.0794415416798357`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.1972245773362196`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.3978952727983707`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.4849066497880004`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.5649493574615367`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.6390573296152584`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.772588722239781, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.833213344056216, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.8903717578961645`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {2.9444389791664403`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {3.6888794541139363`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {4.0943445622221, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}}}, {{{-6.907755278982137, 
       FormBox[
        TagBox[
         InterpretationBox["\"0.001\"", 0.001, AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {-5.298317366548036, 
       FormBox[
        TagBox[
         InterpretationBox["\"0.005\"", 0.005, AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {-4.605170185988092, 
       FormBox[
        TagBox[
         InterpretationBox["\"0.010\"", 0.01, AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {-2.995732273553991, 
       FormBox[
        TagBox[
         InterpretationBox["\"0.050\"", 0.05, AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {-2.302585092994046, 
       FormBox[
        TagBox[
         InterpretationBox["\"0.100\"", 0.1, AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {-0.6931471805599453, 
       FormBox[
        TagBox[
         InterpretationBox["\"0.500\"", 0.5, AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {0., 
       FormBox[
        TagBox[
         InterpretationBox["\"1.000\"", 1., AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {1.6094379124341003`, 
       FormBox[
        TagBox[
         InterpretationBox["\"5.000\"", 5., AutoDelete -> True], 
         NumberForm[#, {
           DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
        TraditionalForm]}, {-6.214608098422191, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-5.809142990314028, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-5.521460917862246, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-5.115995809754082, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-4.961845129926823, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-4.8283137373023015`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-4.710530701645918, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-3.912023005428146, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-3.506557897319982, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-3.2188758248682006`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-2.8134107167600364`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-2.6592600369327783`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-2.5257286443082556`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-2.4079456086518722`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-1.6094379124341003`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-1.203972804325936, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-0.916290731874155, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-0.5108256237659907, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-0.35667494393873245`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-0.2231435513142097, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-0.10536051565782628`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.6931471805599453, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {1.0986122886681098`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {1.3862943611198906`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}}, {{-6.907755278982137, 
       FormBox["\"\"", TraditionalForm]}, {-5.298317366548036, 
       FormBox["\"\"", TraditionalForm]}, {-4.605170185988092, 
       FormBox["\"\"", TraditionalForm]}, {-2.995732273553991, 
       FormBox["\"\"", TraditionalForm]}, {-2.302585092994046, 
       FormBox["\"\"", TraditionalForm]}, {-0.6931471805599453, 
       FormBox["\"\"", TraditionalForm]}, {0., 
       FormBox["\"\"", TraditionalForm]}, {1.6094379124341003`, 
       FormBox["\"\"", TraditionalForm]}, {-6.214608098422191, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-5.809142990314028, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-5.521460917862246, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-5.115995809754082, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-4.961845129926823, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-4.8283137373023015`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-4.710530701645918, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-3.912023005428146, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-3.506557897319982, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-3.2188758248682006`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-2.8134107167600364`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-2.6592600369327783`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-2.5257286443082556`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-2.4079456086518722`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-1.6094379124341003`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-1.203972804325936, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-0.916290731874155, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-0.5108256237659907, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-0.35667494393873245`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-0.2231435513142097, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {-0.10536051565782628`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {0.6931471805599453, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {1.0986122886681098`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}, {1.3862943611198906`, 
       FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
        Thickness[0.001]}}}}},
  GridLines->{None, None},
  Method->{},
  PlotRange->{{-6.907755278982137, 0}, {-0.6931471805599453, 
    0.6931471805599453}},
  PlotRangeClipping->True,
  PlotRangePadding->{
    Scaled[0.02], 
    Scaled[0.02]},
  Ticks->{{{-6.907755278982137, 
      FormBox[
       TagBox[
        InterpretationBox["\"0.001\"", 0.001, AutoDelete -> True], 
        NumberForm[#, {
          DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {-5.298317366548036, 
      FormBox[
       TagBox[
        InterpretationBox["\"0.005\"", 0.005, AutoDelete -> True], 
        NumberForm[#, {
          DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {-4.605170185988092, 
      FormBox[
       TagBox[
        InterpretationBox["\"0.010\"", 0.01, AutoDelete -> True], 
        NumberForm[#, {
          DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {-2.995732273553991, 
      FormBox[
       TagBox[
        InterpretationBox["\"0.050\"", 0.05, AutoDelete -> True], 
        NumberForm[#, {
          DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {-2.302585092994046, 
      FormBox[
       TagBox[
        InterpretationBox["\"0.100\"", 0.1, AutoDelete -> True], 
        NumberForm[#, {
          DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {-0.6931471805599453, 
      FormBox[
       TagBox[
        InterpretationBox["\"0.500\"", 0.5, AutoDelete -> True], 
        NumberForm[#, {
          DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {0., 
      FormBox[
       TagBox[
        InterpretationBox["\"1.000\"", 1., AutoDelete -> True], NumberForm[#, {
          DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {1.6094379124341003`, 
      FormBox[
       TagBox[
        InterpretationBox["\"5.000\"", 5., AutoDelete -> True], NumberForm[#, {
          DirectedInfinity[1], 3}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {-6.214608098422191, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-5.809142990314028, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-5.521460917862246, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-5.115995809754082, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-4.961845129926823, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-4.8283137373023015`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-4.710530701645918, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-3.912023005428146, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-3.506557897319982, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-3.2188758248682006`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-2.8134107167600364`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-2.6592600369327783`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-2.5257286443082556`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-2.4079456086518722`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-1.6094379124341003`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-1.203972804325936, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-0.916290731874155, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-0.5108256237659907, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-0.35667494393873245`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-0.2231435513142097, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {-0.10536051565782628`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {0.6931471805599453, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {1.0986122886681098`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {1.3862943611198906`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}}, {{0., 
      FormBox[
       TagBox[
        InterpretationBox["\"1.0\"", 1., AutoDelete -> True], NumberForm[#, {
          DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {2.302585092994046, 
      FormBox[
       TagBox[
        InterpretationBox["\"10.0\"", 10., AutoDelete -> True], NumberForm[#, {
          DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {1.6094379124341003`, 
      FormBox[
       TagBox[
        InterpretationBox["\"5.0\"", 5., AutoDelete -> True], NumberForm[#, {
          DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {3.912023005428146, 
      FormBox[
       TagBox[
        InterpretationBox["\"50.0\"", 50., AutoDelete -> True], NumberForm[#, {
          DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {0.6931471805599453, 
      FormBox[
       TagBox[
        InterpretationBox["\"2.0\"", 2., AutoDelete -> True], NumberForm[#, {
          DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {2.995732273553991, 
      FormBox[
       TagBox[
        InterpretationBox["\"20.0\"", 20., AutoDelete -> True], NumberForm[#, {
          DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {1.0986122886681098`, 
      FormBox[
       TagBox[
        InterpretationBox["\"3.0\"", 3., AutoDelete -> True], NumberForm[#, {
          DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {3.4011973816621555`, 
      FormBox[
       TagBox[
        InterpretationBox["\"30.0\"", 30., AutoDelete -> True], NumberForm[#, {
          DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {0.4054651081081644, 
      FormBox[
       TagBox[
        InterpretationBox["\"1.5\"", 1.5, AutoDelete -> True], 
        NumberForm[#, {
          DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {2.70805020110221, 
      FormBox[
       TagBox[
        InterpretationBox["\"15.0\"", 15., AutoDelete -> True], NumberForm[#, {
          DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {1.9459101490553132`, 
      FormBox[
       TagBox[
        InterpretationBox["\"7.0\"", 7., AutoDelete -> True], NumberForm[#, {
          DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {4.248495242049359, 
      FormBox[
       TagBox[
        InterpretationBox["\"70.0\"", 70., AutoDelete -> True], NumberForm[#, {
          DirectedInfinity[1], 1}, NumberPadding -> {"", "0"}]& ], 
       TraditionalForm]}, {0.09531017980432493, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {0.1823215567939548, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {0.26236426446749106`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {0.336472236621213, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {0.47000362924573563`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {0.5306282510621705, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {0.5877866649021191, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {0.6418538861723948, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {1.3862943611198906`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {1.791759469228055, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {2.0794415416798357`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {2.1972245773362196`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {2.3978952727983707`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {2.4849066497880004`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {2.5649493574615367`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {2.6390573296152584`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {2.772588722239781, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {2.833213344056216, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {2.8903717578961645`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {2.9444389791664403`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {3.6888794541139363`, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}, {4.0943445622221, 
      FormBox["\"\"", TraditionalForm], {0.00375, 0.}, {
       Thickness[0.001]}}}}]], "Output",
 CellChangeTimes->{{3.592223720814096*^9, 3.5922237762299767`*^9}, 
   3.5922238545593777`*^9, 3.5922239074637327`*^9, 3.592223979082264*^9, 
   3.592772113604527*^9, 3.5927722376036997`*^9, 3.5929147573631077`*^9, 
   3.5931914796729517`*^9, 3.593385029277939*^9, 3.593387546097588*^9, 
   3.59344382090629*^9, 3.593443979026782*^9, 3.593445032682391*^9, 
   3.593446176234935*^9, 3.5934466410575933`*^9, 3.593469401694948*^9, 
   3.5935258156296883`*^9, 3.593525986435121*^9, 3.593963584146945*^9, 
   3.59396466531605*^9, 3.594040334114744*^9, 3.594040853644722*^9, 
   3.5940430469293623`*^9, 3.5944913653943653`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"D", "[", 
  RowBox[{
   RowBox[{"TopHat", "[", 
    RowBox[{"k", ",", " ", "RR"}], "]"}], ",", " ", "RR"}], "]"}]], "Input",
 CellChangeTimes->{{3.592869572639958*^9, 3.5928695813045273`*^9}}],

Cell[BoxData[
 RowBox[{
  SuperscriptBox["TopHat", 
   TagBox[
    RowBox[{"(", 
     RowBox[{"0", ",", "1"}], ")"}],
    Derivative],
   MultilineFunction->None], "[", 
  RowBox[{"k", ",", "RR"}], "]"}]], "Output",
 CellChangeTimes->{3.592869581704977*^9, 3.592914757432687*^9, 
  3.5931914797866697`*^9, 3.5933850293324757`*^9, 3.5933875461475677`*^9, 
  3.593443821359991*^9, 3.5934439797877502`*^9, 3.593445033132846*^9, 
  3.593446176460153*^9, 3.593446641772398*^9, 3.593469401767643*^9, 
  3.5935258171798677`*^9, 3.593525987379024*^9, 3.593963584207533*^9, 
  3.593964665381386*^9, 3.594040336264349*^9, 3.59404085371315*^9, 
  3.594043047046748*^9, 3.594491365464533*^9}]
}, Open  ]]
}, Closed]],

Cell[CellGroupData[{

Cell["Setup integral definitions\[LineSeparator]", "Subsection",
 CellChangeTimes->{{3.59344605007073*^9, 3.59344605264707*^9}, {
   3.5934461558923893`*^9, 3.5934461600687933`*^9}, 3.59396384690987*^9}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"k0", "=", 
    RowBox[{"0.002", "/", "h"}]}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{"In", " ", 
    RowBox[{"h", "/", "Mpc"}]}], " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ns", "=", "0.96"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Lx", "=", "2400"}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{"In", " ", 
    RowBox[{"Mpc", "/", "h"}]}], " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Np", "=", 
   SuperscriptBox["1024", "3"]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Mp", " ", "=", " ", 
    RowBox[{"9.648", " ", 
     SuperscriptBox["10", "11"]}]}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{"mass", " ", "of", " ", "a", " ", "particle", " ", "in", " ", 
    SuperscriptBox["h", 
     RowBox[{"-", "1"}]], 
    SubscriptBox["M", "sun"]}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Rho]bar", " ", "=", " ", 
    FractionBox[
     RowBox[{"Np", "*", "Mp"}], 
     SuperscriptBox[
      RowBox[{"(", "Lx", ")"}], "3"]]}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{"in", " ", 
    SuperscriptBox["h", "2"], 
    RowBox[{
     SubscriptBox["M", "sun"], "/", 
     SuperscriptBox["Mpc", "3"]}]}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"Mass", "[", "R_", "]"}], ":=", " ", 
    RowBox[{
     FractionBox["4", "3"], "\[Pi]", " ", 
     SuperscriptBox["R", "3"], " ", "*", "\[Rho]bar"}]}], ";"}], "   ", 
  RowBox[{"(*", " ", 
   RowBox[{"in", " ", 
    SuperscriptBox["h", 
     RowBox[{"-", "1"}]], 
    SubscriptBox["M", "sun"]}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Aphi", "=", 
   RowBox[{"8.23", " ", 
    SuperscriptBox["10", 
     RowBox[{"-", "10"}]]}]}], ";", 
  RowBox[{"Aphi", "=", 
   SuperscriptBox["10", 
    RowBox[{"-", "5"}]]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"kmin", "=", 
    RowBox[{"2", "*", 
     RowBox[{"\[Pi]", "/", "Lx"}]}]}], ";", " ", 
   RowBox[{"kmax", "=", 
    RowBox[{"1024", "*", "kmin"}]}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{"in", " ", 
    RowBox[{"h", "/", "Mpc"}]}], " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Pvar", "[", "k_", "]"}], ":=", 
   RowBox[{"If", "[", 
    RowBox[{
     RowBox[{"k", ">=", "kmin"}], ",", " ", 
     RowBox[{"Aphi", "*", 
      RowBox[{
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{"k", "/", "k0"}], ")"}], 
        RowBox[{"ns", "-", "1"}]], "/", 
       SuperscriptBox["k", "3"]}]}], ",", " ", "0."}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"TopHat", "[", 
    RowBox[{"k_", ",", "R_"}], "]"}], ":=", 
   FractionBox[
    RowBox[{
     RowBox[{"3", " ", 
      RowBox[{"Sin", "[", 
       RowBox[{"k", " ", "R"}], "]"}]}], "-", 
     RowBox[{"3", 
      RowBox[{"(", 
       RowBox[{"k", " ", "R"}], ")"}], " ", 
      RowBox[{"Cos", "[", 
       RowBox[{"k", " ", "R"}], "]"}]}]}], 
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{"k", " ", "R"}], ")"}], "3"]]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"TopHatD", "[", 
    RowBox[{"k_", ",", "R_"}], "]"}], ":=", 
   RowBox[{
    FractionBox[
     RowBox[{"3", " ", 
      RowBox[{"Sin", "[", 
       RowBox[{"k", " ", "R"}], "]"}]}], 
     RowBox[{"k", " ", 
      SuperscriptBox["R", "2"]}]], "-", 
    FractionBox[
     RowBox[{"3", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"-", "3"}], " ", "k", " ", "R", " ", 
         RowBox[{"Cos", "[", 
          RowBox[{"k", " ", "R"}], "]"}]}], "+", 
        RowBox[{"3", " ", 
         RowBox[{"Sin", "[", 
          RowBox[{"k", " ", "R"}], "]"}]}]}], ")"}]}], 
     RowBox[{
      SuperscriptBox["k", "3"], " ", 
      SuperscriptBox["R", "4"]}]]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Alpha]", "[", "k_", "]"}], ":=", 
   RowBox[{"MEH", "[", 
    RowBox[{"k", "*", "h"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"\[Alpha]M", "[", 
     RowBox[{"k_", ",", "R_"}], "]"}], ":=", 
    RowBox[{
     RowBox[{"MEH", "[", 
      RowBox[{"k", "*", "h"}], "]"}], "*", 
     RowBox[{"TopHat", "[", 
      RowBox[{"k", ",", "R"}], "]"}]}]}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"MEH", "[", "k", "]"}], " ", "takes", " ", "k", " ", "in", " ", 
    SuperscriptBox["Mpc", 
     RowBox[{"-", "1"}]]}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"CosTheta", "[", 
    RowBox[{"\[Phi]1_", ",", "\[Phi]2_", ",", "\[Mu]1_", ",", "\[Mu]2_"}], 
    "]"}], ":=", 
   RowBox[{
    RowBox[{"\[Mu]1", "*", "\[Mu]2"}], "+", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"Cos", "[", "\[Phi]1", "]"}], "*", 
        RowBox[{"Cos", "[", "\[Phi]2", "]"}]}], "+", 
       RowBox[{
        RowBox[{"Sin", "[", "\[Phi]1", "]"}], "*", 
        RowBox[{"Sin", "[", "\[Phi]2", "]"}]}]}], ")"}], "*", 
     RowBox[{"Sqrt", "[", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "-", 
         RowBox[{"\[Mu]1", "^", "2"}]}], ")"}], "*", 
       RowBox[{"(", 
        RowBox[{"1", "-", 
         RowBox[{"\[Mu]2", "^", "2"}]}], ")"}]}], "]"}]}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"I32Value", "[", 
    RowBox[{
    "k_", ",", "R_", ",", " ", "q2_", ",", " ", "\[Mu]2_", ",", " ", 
     "\[Phi]2_", ",", " ", "q1_", ",", " ", "\[Mu]1_", ",", " ", "\[Phi]1_"}],
     "]"}], ":=", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       SuperscriptBox["q2", "2"], "*", 
       RowBox[{"\[Alpha]M", "[", 
        RowBox[{"q2", ",", "R"}], "]"}], "*", 
       RowBox[{"\[Alpha]M", "[", 
        RowBox[{
         RowBox[{"Sqrt", "[", 
          RowBox[{
           SuperscriptBox["q2", "2"], "+", 
           SuperscriptBox["k", "2"], "+", 
           RowBox[{"2", "*", "q2", "*", "k", "*", "\[Mu]2"}]}], "]"}], ",", 
         "R"}], "]"}], "*", 
       SuperscriptBox["q1", "2"], "*", 
       RowBox[{"Pvar", "[", "q1", "]"}], 
       RowBox[{"Pvar", "[", 
        RowBox[{"Sqrt", "[", 
         RowBox[{
          SuperscriptBox["q2", "2"], "+", 
          SuperscriptBox["q1", "2"], "-", 
          RowBox[{"2", "*", "q2", "*", "q1", "*", 
           RowBox[{"CosTheta", "[", 
            RowBox[{
            "\[Phi]1", ",", " ", "\[Phi]2", ",", " ", "\[Mu]1", ",", " ", 
             "\[Mu]2"}], "]"}]}]}], "]"}], "]"}], "*", 
       RowBox[{"(", 
        RowBox[{"Pvar", "[", 
         RowBox[{"Sqrt", "[", 
          RowBox[{
           SuperscriptBox["k", "2"], "+", 
           SuperscriptBox["q1", "2"], "+", 
           RowBox[{"2", "*", "k", "*", "q1", "*", "\[Mu]1"}]}], "]"}], "]"}], 
        ")"}]}], ")"}], "/", "8."}], "/", 
    SuperscriptBox["\[Pi]", "2"]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"I31Value", "[", 
    RowBox[{"k_", ",", "R_", ",", "q1_", ",", "\[Mu]1_"}], "]"}], ":=", 
   RowBox[{
    SuperscriptBox["q1", "2"], "*", 
    RowBox[{"\[Alpha]M", "[", 
     RowBox[{"q1", ",", "R"}], "]"}], "*", 
    RowBox[{"\[Alpha]M", "[", 
     RowBox[{
      RowBox[{"Sqrt", "[", 
       RowBox[{
        SuperscriptBox["q1", "2"], "+", 
        SuperscriptBox["k", "2"], "+", 
        RowBox[{"2", "*", "q1", "*", "k", "*", "\[Mu]1"}]}], "]"}], ",", 
      "R"}], "]"}], "*", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"Pvar", "[", "q1", "]"}], "*", 
         RowBox[{"Pvar", "[", 
          RowBox[{"Sqrt", "[", 
           RowBox[{
            SuperscriptBox["q1", "2"], "+", 
            SuperscriptBox["k", "2"], "+", 
            RowBox[{"2", "*", "q1", "*", "k", "*", "\[Mu]1"}]}], "]"}], 
          "]"}]}], ")"}], "+", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"Pvar", "[", 
          RowBox[{"Sqrt", "[", 
           RowBox[{
            SuperscriptBox["q1", "2"], "+", 
            SuperscriptBox["k", "2"], "+", 
            RowBox[{"2", "*", "q1", "*", "k", "*", "\[Mu]1"}]}], "]"}], "]"}],
          "*", 
         RowBox[{"Pvar", "[", "k", "]"}]}], ")"}], "+", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"Pvar", "[", "k", "]"}], "*", 
         RowBox[{"Pvar", "[", "q1", "]"}]}], ")"}]}], ")"}], "/", "4.0"}]}]}],
   ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"ISigSq", "[", 
     RowBox[{"k_", ",", "R_"}], "]"}], ":=", 
    RowBox[{
     SuperscriptBox[
      RowBox[{"\[Alpha]M", "[", 
       RowBox[{"k", ",", "R"}], "]"}], "2"], "*", 
     RowBox[{"Pvar", "[", "k", "]"}], "*", "k", "*", "k"}]}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{
    "for", " ", "SigmaMSq", " ", "computation", " ", "to", " ", "get", " ", 
     "A", " ", "such", " ", "that", " ", 
     SubscriptBox["\[Sigma]", "8"]}], "=", "0.8"}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     RowBox[{",", " ", 
      RowBox[{"{", 
       RowBox[{"q2", ",", " ", "kmin", ",", " ", "kmax"}], "}"}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"\[Mu]2", ",", " ", 
        RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"\[Phi]2", ",", " ", "0", ",", " ", 
        RowBox[{"2", "*", "\[Pi]"}]}], "}"}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"q1", ",", " ", "kmin", ",", " ", "kmax"}], "}"}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"\[Mu]1", ",", " ", 
        RowBox[{"-", "1"}], ",", " ", "1"}], "}"}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"\[Phi]1", ",", " ", "0", ",", " ", 
        RowBox[{"2", "*", "\[Pi]"}]}], "}"}]}], "]"}], ";"}], "*)"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"I31DLog", "[", 
     RowBox[{"k_", ",", "R_", ",", "q1_", ",", "\[Mu]1_"}], "]"}], ":=", 
    RowBox[{
     SuperscriptBox["q1", "2"], "*", 
     RowBox[{"\[Alpha]", "[", "q1", "]"}], "*", 
     RowBox[{"\[Alpha]", "[", 
      RowBox[{"Sqrt", "[", 
       RowBox[{
        SuperscriptBox["q1", "2"], "+", 
        SuperscriptBox["k", "2"], "+", 
        RowBox[{"2", "*", "q1", "*", "k", "*", "\[Mu]1"}]}], "]"}], "]"}], 
     "*", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"TopHatD", "[", 
         RowBox[{"q1", ",", "R"}], "]"}], "*", 
        RowBox[{"TopHat", "[", 
         RowBox[{
          RowBox[{"Sqrt", "[", 
           RowBox[{
            SuperscriptBox["q1", "2"], "+", 
            SuperscriptBox["k", "2"], "+", 
            RowBox[{"2", "*", "q1", "*", "k", "*", "\[Mu]1"}]}], "]"}], ",", 
          "R"}], "]"}]}], "+", 
       RowBox[{
        RowBox[{"TopHat", "[", 
         RowBox[{"q1", ",", "R"}], "]"}], "*", 
        RowBox[{"TopHatD", "[", 
         RowBox[{
          RowBox[{"Sqrt", "[", 
           RowBox[{
            SuperscriptBox["q1", "2"], "+", 
            SuperscriptBox["k", "2"], "+", 
            RowBox[{"2", "*", "q1", "*", "k", "*", "\[Mu]1"}]}], "]"}], ",", 
          "R"}], "]"}]}]}], ")"}], "*", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Pvar", "[", "q1", "]"}], "*", 
          RowBox[{"Pvar", "[", 
           RowBox[{"Sqrt", "[", 
            RowBox[{
             SuperscriptBox["q1", "2"], "+", 
             SuperscriptBox["k", "2"], "+", 
             RowBox[{"2", "*", "q1", "*", "k", "*", "\[Mu]1"}]}], "]"}], 
           "]"}]}], ")"}], "+", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Pvar", "[", 
           RowBox[{"Sqrt", "[", 
            RowBox[{
             SuperscriptBox["q1", "2"], "+", 
             SuperscriptBox["k", "2"], "+", 
             RowBox[{"2", "*", "q1", "*", "k", "*", "\[Mu]1"}]}], "]"}], 
           "]"}], "*", 
          RowBox[{"Pvar", "[", "k", "]"}]}], ")"}], "+", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Pvar", "[", "k", "]"}], "*", 
          RowBox[{"Pvar", "[", "q1", "]"}]}], ")"}]}], ")"}], "/", 
      "4.0"}]}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"I32DLog", "[", 
    RowBox[{
    "k_", ",", "R_", ",", " ", "q2_", ",", " ", "\[Mu]2_", ",", " ", 
     "\[Phi]2_", ",", " ", "q1_", ",", " ", "\[Mu]1_", ",", " ", "\[Phi]1_"}],
     "]"}], ":=", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       SuperscriptBox["q2", "2"], "*", 
       RowBox[{"\[Alpha]", "[", "q2", "]"}], "*", 
       RowBox[{"\[Alpha]", "[", 
        RowBox[{"Sqrt", "[", 
         RowBox[{
          SuperscriptBox["q2", "2"], "+", 
          SuperscriptBox["k", "2"], "+", 
          RowBox[{"2", "*", "q2", "*", "k", "*", "\[Mu]2"}]}], "]"}], "]"}], 
       "*", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{"TopHatD", "[", 
           RowBox[{"q2", ",", "R"}], "]"}], "*", 
          RowBox[{"TopHat", "[", 
           RowBox[{
            RowBox[{"Sqrt", "[", 
             RowBox[{
              SuperscriptBox["q2", "2"], "+", 
              SuperscriptBox["k", "2"], "+", 
              RowBox[{"2", "*", "q2", "*", "k", "*", "\[Mu]2"}]}], "]"}], ",",
             "R"}], "]"}]}], "+", 
         RowBox[{
          RowBox[{"TopHat", "[", 
           RowBox[{"q2", ",", "R"}], "]"}], "*", 
          RowBox[{"TopHatD", "[", 
           RowBox[{
            RowBox[{"Sqrt", "[", 
             RowBox[{
              SuperscriptBox["q2", "2"], "+", 
              SuperscriptBox["k", "2"], "+", 
              RowBox[{"2", "*", "q2", "*", "k", "*", "\[Mu]2"}]}], "]"}], ",",
             "R"}], "]"}]}]}], ")"}], "*", 
       SuperscriptBox["q1", "2"], "*", 
       RowBox[{"Pvar", "[", "q1", "]"}], 
       RowBox[{"Pvar", "[", 
        RowBox[{"Sqrt", "[", 
         RowBox[{
          SuperscriptBox["q2", "2"], "+", 
          SuperscriptBox["q1", "2"], "-", 
          RowBox[{"2", "*", "q2", "*", "q1", "*", 
           RowBox[{"CosTheta", "[", 
            RowBox[{
            "\[Phi]1", ",", " ", "\[Phi]2", ",", " ", "\[Mu]1", ",", " ", 
             "\[Mu]2"}], "]"}]}]}], "]"}], "]"}], "*", 
       RowBox[{"(", 
        RowBox[{"Pvar", "[", 
         RowBox[{"Sqrt", "[", 
          RowBox[{
           SuperscriptBox["k", "2"], "+", 
           SuperscriptBox["q1", "2"], "+", 
           RowBox[{"2", "*", "k", "*", "q1", "*", "\[Mu]1"}]}], "]"}], "]"}], 
        ")"}]}], ")"}], "/", "8."}], "/", 
    SuperscriptBox["\[Pi]", "2"]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"I42Value", "[", 
    RowBox[{
    "k_", ",", "R_", ",", " ", "q2_", ",", " ", "\[Mu]2_", ",", " ", 
     "\[Phi]2_", ",", " ", "q1_", ",", " ", "\[Mu]1_", ",", " ", "\[Phi]1_", 
     ",", " ", "p_", ",", " ", "\[Mu]_", ",", " ", "\[Phi]_"}], "]"}], ":=", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       SuperscriptBox["q1", "2"], "*", 
       SuperscriptBox["q2", "2"], "*", 
       SuperscriptBox["p", "2"], "*", 
       SuperscriptBox[
        RowBox[{"\[Alpha]M", "[", 
         RowBox[{"q1", ",", "R"}], "]"}], "2"], "*", 
       SuperscriptBox[
        RowBox[{"\[Alpha]M", "[", 
         RowBox[{"q2", ",", "R"}], "]"}], "2"], "*", 
       RowBox[{"Pvar", "[", "p", "]"}], "*", 
       RowBox[{"Pvar", "[", 
        RowBox[{"Sqrt", "[", 
         RowBox[{
          SuperscriptBox["q1", "2"], "-", 
          RowBox[{"2", "*", "q1", "*", "p", "*", 
           RowBox[{"CosTheta", "[", 
            RowBox[{"\[Phi]1", ",", "\[Phi]", ",", "\[Mu]1", ",", "\[Mu]"}], 
            "]"}]}], "+", 
          SuperscriptBox["p", "2"]}], "]"}], "]"}], "*", 
       RowBox[{"Pvar", "[", 
        RowBox[{"Sqrt", "[", 
         RowBox[{
          SuperscriptBox["k", "2"], "+", 
          SuperscriptBox["q1", "2"], "+", 
          SuperscriptBox["p", "2"], "-", 
          RowBox[{"2", "*", "k", "*", "q1", "*", "\[Mu]1"}], "+", 
          RowBox[{"2", "*", "k", "*", "p", "*", "\[Mu]"}], "-", 
          RowBox[{"2", "*", "q1", "*", "p", "*", 
           RowBox[{"CosTheta", "[", 
            RowBox[{"\[Phi]1", ",", "\[Phi]", ",", "\[Mu]1", ",", "\[Mu]"}], 
            "]"}]}]}], "]"}], "]"}], "*", 
       RowBox[{"Pvar", "[", 
        RowBox[{"Sqrt", "[", 
         RowBox[{
          SuperscriptBox["p", "2"], "+", 
          SuperscriptBox["q1", "2"], "+", 
          SuperscriptBox["q2", "2"], "-", 
          RowBox[{"2", "*", "p", "*", "q1", "*", 
           RowBox[{"CosTheta", "[", 
            RowBox[{"\[Phi]", ",", "\[Phi]1", ",", "\[Mu]", ",", "\[Mu]1"}], 
            "]"}]}], "-", 
          RowBox[{"2", "*", "p", "*", "q2", "*", 
           RowBox[{"CosTheta", "[", 
            RowBox[{"\[Phi]", ",", "\[Phi]2", ",", "\[Mu]", ",", "\[Mu]2"}], 
            "]"}]}], "+", 
          RowBox[{"2", "*", "q1", "*", "q2", "*", 
           RowBox[{"CosTheta", "[", 
            RowBox[{"\[Phi]1", ",", "\[Phi]2", ",", "\[Mu]1", ",", "\[Mu]2"}],
             "]"}]}]}], "]"}], "]"}]}], ")"}], "/", "64."}], "/", 
    SuperscriptBox["\[Pi]", "3"]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"I42Approx", "[", 
    RowBox[{
    "k_", ",", "R_", ",", "p1_", ",", " ", "\[Mu]1_", ",", " ", "\[Phi]1_", 
     ",", " ", "p_", ",", " ", "\[Mu]_", ",", " ", "\[Phi]_"}], "]"}], ":=", 
   RowBox[{
    SuperscriptBox["p1", "2"], "*", 
    SuperscriptBox["p", "2"], "*", 
    SuperscriptBox[
     RowBox[{"\[Alpha]M", "[", 
      RowBox[{"p1", ",", "R"}], "]"}], "2"], "*", 
    RowBox[{"Pvar", "[", "p", "]"}], "*", 
    RowBox[{"Pvar", "[", 
     RowBox[{"Sqrt", "[", 
      RowBox[{
       SuperscriptBox["p1", "2"], "+", 
       SuperscriptBox["p", "2"], "-", 
       RowBox[{"2", "*", "p1", "*", "p", "*", 
        RowBox[{"CosTheta", "[", 
         RowBox[{"\[Phi]1", ",", "\[Phi]", ",", "\[Mu]1", ",", "\[Mu]"}], 
         "]"}]}]}], "]"}], "]"}], "*", 
    RowBox[{
     RowBox[{
      RowBox[{"Pvar", "[", 
       RowBox[{"Sqrt", "[", 
        RowBox[{
         SuperscriptBox["k", "2"], "+", 
         SuperscriptBox["p1", "2"], "+", 
         SuperscriptBox["p", "2"], "-", 
         RowBox[{"2", "*", "k", "*", "p1", "*", "\[Mu]1"}], "+", 
         RowBox[{"2", "*", "k", "*", "p", "*", "\[Mu]"}], "-", 
         RowBox[{"2", "*", "p1", "*", "p", "*", 
          RowBox[{"CosTheta", "[", 
           RowBox[{
           "\[Phi]1", ",", "\[Phi]", ",", "\[Mu]1", ",", " ", "\[Mu]"}], 
           "]"}]}]}], "]"}], "]"}], "/", "16."}], "/", 
     SuperscriptBox["\[Pi]", "2"]}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"B4Exact", "[", 
    RowBox[{
    "k_", ",", " ", "p1_", ",", " ", "p1\[Mu]_", ",", "p1\[Phi]_", ",", "p2_",
      ",", "p2\[Mu]_", ",", "p2\[Phi]_", ",", "p_", ",", "\[Mu]_", ",", 
     "\[Phi]_"}], "]"}], ":=", 
   RowBox[{
    SuperscriptBox["p", "2"], 
    RowBox[{"Pvar", "[", "p", "]"}], "*", 
    RowBox[{"Pvar", "[", 
     RowBox[{"Sqrt", "[", 
      RowBox[{
       SuperscriptBox["p1", "2"], "+", 
       SuperscriptBox["p", "2"], "-", 
       RowBox[{"2", "*", "p1", "*", "p", "*", 
        RowBox[{"CosTheta", "[", 
         RowBox[{"p1\[Phi]", ",", "\[Phi]", ",", "p1\[Mu]", ",", "\[Mu]"}], 
         "]"}]}]}], "]"}], "]"}], "*", 
    RowBox[{"Pvar", "[", 
     RowBox[{"Sqrt", "[", 
      RowBox[{
       SuperscriptBox["k", "2"], "+", 
       SuperscriptBox["p1", "2"], "+", 
       SuperscriptBox["p", "2"], "-", 
       RowBox[{"2", "*", "k", "*", "p1", "*", "p1\[Mu]"}], "+", 
       RowBox[{"2", "*", "k", "*", "p", "*", "\[Mu]"}], "-", 
       RowBox[{"2", "*", "p1", "*", "p", "*", 
        RowBox[{"CosTheta", "[", 
         RowBox[{"p1\[Phi]", ",", "\[Phi]", ",", "p1\[Mu]", ",", "\[Mu]"}], 
         "]"}]}]}], "]"}], "]"}], "*", 
    RowBox[{"Pvar", "[", 
     RowBox[{"Sqrt", "[", 
      RowBox[{
       SuperscriptBox["p", "2"], "+", 
       SuperscriptBox["p1", "2"], "+", 
       SuperscriptBox["p2", "2"], "-", 
       RowBox[{"2", "*", "p", "*", "p1", "*", 
        RowBox[{"CosTheta", "[", 
         RowBox[{"\[Phi]", ",", "p1\[Phi]", ",", "\[Mu]", ",", "p1\[Mu]"}], 
         "]"}]}], "-", 
       RowBox[{"2", "*", "p", "*", "p2", "*", 
        RowBox[{"CosTheta", "[", 
         RowBox[{"\[Phi]", ",", "p2\[Phi]", ",", "\[Mu]", ",", "p2\[Mu]"}], 
         "]"}]}], "+", 
       RowBox[{"2", "*", "p1", "*", "p2", "*", 
        RowBox[{"CosTheta", "[", 
         RowBox[{
         "p1\[Phi]", ",", "p2\[Phi]", ",", "p1\[Mu]", ",", "p2\[Mu]"}], 
         "]"}]}]}], "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"B4Approx", "[", 
   RowBox[{
   "k_", ",", " ", "p1_", ",", " ", "p1\[Mu]_", ",", "p1\[Phi]_", ",", "p2_", 
    ",", "p_", ",", "\[Mu]_", ",", "\[Phi]_"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"Pvar", "[", "p2", "]"}], "*", 
   SuperscriptBox["p", "2"], 
   RowBox[{"Pvar", "[", "p", "]"}], "*", 
   RowBox[{"Pvar", "[", 
    RowBox[{"Sqrt", "[", 
     RowBox[{
      SuperscriptBox["p1", "2"], "+", 
      SuperscriptBox["p", "2"], "-", 
      RowBox[{"2", "*", "p1", "*", "p", "*", 
       RowBox[{"CosTheta", "[", 
        RowBox[{"p1\[Phi]", ",", "\[Phi]", ",", "p1\[Mu]", ",", "\[Mu]"}], 
        "]"}]}]}], "]"}], "]"}], "*", 
   RowBox[{"Pvar", "[", 
    RowBox[{"Sqrt", "[", 
     RowBox[{
      SuperscriptBox["k", "2"], "+", 
      SuperscriptBox["p1", "2"], "+", 
      SuperscriptBox["p", "2"], "-", 
      RowBox[{"2", "*", "k", "*", "p1", "*", "p1\[Mu]"}], "+", 
      RowBox[{"2", "*", "k", "*", "p", "*", "\[Mu]"}], "-", 
      RowBox[{"2", "*", "p1", "*", "p", "*", 
       RowBox[{"CosTheta", "[", 
        RowBox[{"p1\[Phi]", ",", "\[Phi]", ",", "p1\[Mu]", ",", "\[Mu]"}], 
        "]"}]}]}], "]"}], "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"IPPApprox", "[", 
    RowBox[{"k_", ",", " ", "p_", ",", "\[Mu]_"}], "]"}], ":=", 
   RowBox[{
    SuperscriptBox["p", "2"], "*", 
    RowBox[{"Pvar", "[", "p", "]"}], "*", 
    RowBox[{"Pvar", "[", 
     RowBox[{"Sqrt", "[", 
      RowBox[{
       SuperscriptBox["p", "2"], "+", 
       SuperscriptBox["k", "2"], "+", 
       RowBox[{"2", "*", "k", "*", "p", "*", "\[Mu]"}]}], "]"}], "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"IPPInt", "[", "k_", "]"}], ":=", 
   RowBox[{"Suave", "[", 
    RowBox[{
     RowBox[{"IPPApprox", "[", 
      RowBox[{"k", ",", "p", ",", "\[Mu]"}], "]"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"p", ",", " ", "kmin", ",", " ", "kmax"}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"\[Mu]", ",", " ", 
       RowBox[{"-", "1."}], ",", " ", "1"}], "}"}], ",", " ", 
     RowBox[{"Verbose", "\[Rule]", "0"}], ",", " ", 
     RowBox[{"MaxPoints", "\[Rule]", "1000000"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"SigmapsiSq", "[", "R_", "]"}], ":=", 
   RowBox[{"Suave", "[", 
    RowBox[{
     RowBox[{"ISigSq", "[", 
      RowBox[{"k", ",", " ", "R"}], "]"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"k", ",", " ", "kmin", ",", " ", "kmax"}], "}"}], ",", " ", 
     RowBox[{"Verbose", "\[Rule]", "0"}], ",", " ", 
     RowBox[{"MaxPoints", "\[Rule]", "100000"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.5920663919167843`*^9, 3.592066396532892*^9}, {
   3.592066458555657*^9, 3.592066805894929*^9}, 3.592067410093257*^9, {
   3.592067448640141*^9, 3.592067753453969*^9}, {3.592067798100644*^9, 
   3.5920679169902067`*^9}, {3.592068001218309*^9, 3.592068040062264*^9}, 
   3.592069337368183*^9, {3.59206945968678*^9, 3.592069460103941*^9}, {
   3.5920695056602287`*^9, 3.5920695595299063`*^9}, {3.592069608487071*^9, 
   3.592069751026957*^9}, {3.592150839910246*^9, 3.59215084780333*^9}, {
   3.592153537485997*^9, 3.592153566789893*^9}, {3.592223259730954*^9, 
   3.592223283894908*^9}, {3.5922255267657967`*^9, 3.5922255284008102`*^9}, {
   3.5927612547253933`*^9, 3.592761335605227*^9}, {3.592761494875949*^9, 
   3.592761500818406*^9}, {3.592762052046629*^9, 3.592762121015386*^9}, {
   3.5927622925193157`*^9, 3.592762294594532*^9}, {3.592763655449903*^9, 
   3.5927638177349157`*^9}, {3.592764111914147*^9, 3.592764143295429*^9}, {
   3.5927710299524918`*^9, 3.592771045323942*^9}, {3.592821990166234*^9, 
   3.592821999012002*^9}, {3.59282217043034*^9, 3.592822174090159*^9}, 
   3.59282628851263*^9, {3.592827202419595*^9, 3.5928272029499607`*^9}, {
   3.592827556201213*^9, 3.5928275568993597`*^9}, {3.59282762194914*^9, 
   3.592827650945356*^9}, {3.592850891018014*^9, 3.592850918549065*^9}, {
   3.592851016030954*^9, 3.5928510275058603`*^9}, {3.592851057811228*^9, 
   3.592851118888485*^9}, 3.5928514483024893`*^9, {3.592869530552868*^9, 
   3.592869554335053*^9}, {3.5928695886175547`*^9, 3.592869708033636*^9}, {
   3.592871608495091*^9, 3.5928716798080177`*^9}, {3.593191589501388*^9, 
   3.593191635049176*^9}, {3.593191724405734*^9, 3.593191742744581*^9}, {
   3.593191837638856*^9, 3.5931919100235767`*^9}, {3.593192155329774*^9, 
   3.593192156225503*^9}, {3.593192232513816*^9, 3.593192262598289*^9}, {
   3.593192295004855*^9, 3.593192419986742*^9}, {3.593192490670388*^9, 
   3.593192493915887*^9}, {3.593192534582417*^9, 3.5931925356987762`*^9}, {
   3.59319856535787*^9, 3.593198582837366*^9}, {3.5932197161334267`*^9, 
   3.5932197236399603`*^9}, 3.593276882937324*^9, {3.593387291530457*^9, 
   3.5933872998572598`*^9}, {3.593387330264431*^9, 3.593387340145187*^9}, {
   3.5934710587910433`*^9, 3.593471061909923*^9}, {3.5935158998259363`*^9, 
   3.593515910618945*^9}, {3.593516143028067*^9, 3.593516143703466*^9}, {
   3.593964025405846*^9, 3.593964075055188*^9}, {3.5939641880151978`*^9, 
   3.593964590945448*^9}, {3.5940488369145308`*^9, 3.594048878431732*^9}, {
   3.594049457664654*^9, 3.594049515793841*^9}, {3.5940495734028273`*^9, 
   3.594049715478744*^9}, {3.5940500291888514`*^9, 3.594050063497674*^9}, {
   3.594050257334449*^9, 3.594050350324565*^9}, {3.594491421811514*^9, 
   3.594491445469487*^9}, {3.594491506340189*^9, 3.5944915392165117`*^9}, {
   3.59449162133358*^9, 3.594491728224533*^9}, {3.5944918403478394`*^9, 
   3.594491904185679*^9}, {3.594492289397223*^9, 3.594492291995483*^9}, {
   3.594501823743679*^9, 3.594501826683626*^9}, {3.594502876084482*^9, 
   3.594502877502033*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"klist", " ", "=", " ", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"N", "[", "kmin", "]"}], "+", 
       RowBox[{
        FractionBox[
         RowBox[{"N", "[", "kmin", "]"}], "2"], "*", 
        RowBox[{"(", 
         RowBox[{"i", "-", "1"}], ")"}]}]}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"i", ",", " ", "1", ",", " ", "NRS"}], "}"}]}], "]"}]}], ";"}],
   " ", "//", "N"}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"F32Approx", "=", 
   RowBox[{"ParallelTable", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"klist", "[", 
        RowBox[{"[", "i", "]"}], "]"}], ",", " ", 
       RowBox[{"IPPInt", "[", 
        RowBox[{"klist", "[", 
         RowBox[{"[", "i", "]"}], "]"}], "]"}]}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"i", ",", " ", "1", ",", " ", 
       RowBox[{"Length", "[", "klist", "]"}]}], "}"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.594502050984869*^9, 3.5945020905581636`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"F32List", "=", " ", 
   RowBox[{"ParallelTable", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"klist", "[", 
        RowBox[{"[", "i", "]"}], "]"}], ",", " ", 
       RowBox[{
        RowBox[{
         RowBox[{"IntegratedValue32", "[", 
          RowBox[{
           RowBox[{"klist", "[", 
            RowBox[{"[", "i", "]"}], "]"}], ",", " ", "8."}], "]"}], "[", 
         RowBox[{"[", "1", "]"}], "]"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"i", ",", " ", "1", ",", " ", 
       RowBox[{"Length", "[", "klist", "]"}]}], "}"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.5945018284496613`*^9, 3.594501834830152*^9}, {
  3.594501928572199*^9, 3.594501963103788*^9}, {3.5945020189467916`*^9, 
  3.594502020835885*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(531295\\) function evaluations on \
\\!\\(531\\) subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.594502535880496*^9}],

Cell[BoxData["$Aborted"], "Output",
 CellChangeTimes->{3.594502017026935*^9, 3.594502664348465*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"0.000486627", "*", "7776.23"}]], "Input",
 CellChangeTimes->{{3.594492200970111*^9, 3.594492225363246*^9}, {
  3.5944923046529503`*^9, 3.5944923129023438`*^9}}],

Cell[BoxData["3.78412347621`"], "Output",
 CellChangeTimes->{3.594492225747652*^9, 3.594492313273138*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"IPPInt", "[", "0.01", "]"}]], "Input",
 CellChangeTimes->{{3.594492186822632*^9, 3.5944921918941298`*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \[NoBreak]\\!\\(77038\\)\[NoBreak] function \
evaluations on \[NoBreak]\\!\\(77\\)\[NoBreak] subregions.\"\>"}]], "Message",\
 "MSG",
 CellChangeTimes->{3.5944921953340273`*^9, 3.5944922996200542`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"{", 
   RowBox[{
   "0.00048662655738277116`", ",", "4.79178873998302`*^-7", ",", "1.`"}], 
   "}"}], "}"}]], "Output",
 CellChangeTimes->{3.594492195335705*^9, 3.59449229962155*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"SigmapsiSq", "[", "8.", "]"}]], "Input",
 CellChangeTimes->{{3.5944917776898003`*^9, 3.594491785759767*^9}, {
  3.594491889261241*^9, 3.594491893137168*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \[NoBreak]\\!\\(10000\\)\[NoBreak] function \
evaluations on \[NoBreak]\\!\\(10\\)\[NoBreak] subregions.\"\>"}]], "Message",\
 "MSG",
 CellChangeTimes->{{3.594491894510284*^9, 3.5944919067508907`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"{", 
   RowBox[{
   "7776.228260730612`", ",", "7.368595046230659`", ",", 
    "0.5771531696643812`"}], "}"}], "}"}]], "Output",
 CellChangeTimes->{{3.59449178165121*^9, 3.594491786424679*^9}, {
  3.594491894512426*^9, 3.594491906753355*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"IntegratedValue32", "[", 
  RowBox[{"0.01", ",", " ", "8."}], "]"}]], "Input",
 CellChangeTimes->{{3.594491705162348*^9, 3.5944917535272293`*^9}, {
  3.594491925324304*^9, 3.594491958781393*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \[NoBreak]\\!\\(548213\\)\[NoBreak] function \
evaluations on \[NoBreak]\\!\\(548\\)\[NoBreak] subregions.\"\>"}]], "Message",\
 "MSG",
 CellChangeTimes->{{3.594491709953676*^9, 3.594491757857732*^9}, 
   3.594492156077175*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"{", 
   RowBox[{"3.8100931539380234`", ",", "0.0038089520746813384`", ",", "1.`"}],
    "}"}], "}"}]], "Output",
 CellChangeTimes->{{3.5944917099546747`*^9, 3.5944917578604183`*^9}, 
   3.594492156079401*^9}]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"B4E", "[", 
    RowBox[{
    "k_", ",", "p1_", ",", "p1\[Mu]_", ",", "p1\[Phi]_", ",", "p2_", ",", 
     "p2\[Mu]_", ",", "p2\[Phi]_"}], "]"}], ":=", " ", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"CosTheta", "[", 
      RowBox[{
      "p1\[Phi]", ",", " ", "p2\[Phi]", ",", " ", "p1\[Mu]", ",", " ", 
       "p2\[Mu]"}], "]"}], ",", 
     RowBox[{
      RowBox[{
       RowBox[{"Suave", "[", 
        RowBox[{
         RowBox[{"B4Exact", "[", 
          RowBox[{
          "k", ",", "p1", ",", "p1\[Mu]", ",", "p1\[Phi]", ",", "p2", ",", 
           "p2\[Mu]", ",", "p2\[Phi]", ",", "p", ",", "\[Mu]", ",", 
           "\[Phi]"}], "]"}], ",", " ", 
         RowBox[{"{", 
          RowBox[{"p", ",", " ", "kmin", ",", " ", "kmax"}], "}"}], ",", " ", 
         
         RowBox[{"{", 
          RowBox[{"\[Mu]", ",", 
           RowBox[{"-", "1."}], ",", "1."}], "}"}], ",", " ", 
         RowBox[{"{", 
          RowBox[{"\[Phi]", ",", "0.", ",", 
           RowBox[{"2.", "*", "\[Pi]"}]}], "}"}], ",", " ", 
         RowBox[{"Verbose", "\[Rule]", "0"}], ",", " ", 
         RowBox[{"MaxPoints", "\[Rule]", "1000000"}]}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}], "[", 
      RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"B4A", "[", 
    RowBox[{"k_", ",", "p1_", ",", "p1\[Mu]_", ",", "p1\[Phi]_", ",", "p2_"}],
     "]"}], ":=", " ", 
   RowBox[{"Suave", "[", 
    RowBox[{
     RowBox[{"B4Approx", "[", 
      RowBox[{
      "k", ",", "p1", ",", "p1\[Mu]", ",", "p1\[Phi]", ",", "p2", ",", "p", 
       ",", "\[Mu]", ",", "\[Phi]"}], "]"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"p", ",", " ", "kmin", ",", " ", "kmax"}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"\[Mu]", ",", 
       RowBox[{"-", "1."}], ",", "1."}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"\[Phi]", ",", "0.", ",", 
       RowBox[{"2.", "*", "\[Pi]"}]}], "}"}], ",", " ", 
     RowBox[{"Verbose", "\[Rule]", "0"}], ",", " ", 
     RowBox[{"MaxPoints", "\[Rule]", "1000000"}]}], "]"}]}], ";"}]}], "Input",\

 CellChangeTimes->{{3.593964693631995*^9, 3.593964877023055*^9}, {
   3.59396495924536*^9, 3.5939649629965773`*^9}, {3.594041338356142*^9, 
   3.5940413434681892`*^9}, {3.594041600466813*^9, 3.5940416725921917`*^9}, 
   3.594041744938189*^9, 3.5940420999115*^9, {3.594042436489994*^9, 
   3.594042437962332*^9}, {3.594042966123971*^9, 3.594042969157795*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{
    "lets", " ", "work", " ", "out", " ", "this", " ", "exact", " ", "vs", 
     " ", "approx", " ", "for", " ", "a", " ", "few", " ", "cases"}], ",", 
    " ", 
    RowBox[{
     RowBox[{"say", " ", "k"}], "=", 
     RowBox[{"0.01", " ", "and", " ", "p1"}]}], ",", " ", 
    RowBox[{
     RowBox[{"p2", " ", "~", " ", "0.1"}], ";", " ", 
     RowBox[{
     "then", " ", "we", " ", "have", " ", "to", " ", "choose", " ", "some", 
      " ", "orientation", " ", "of", " ", 
      RowBox[{"p1", "."}]}]}]}], " ", "*)"}], 
  "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.593964922107875*^9, 3.593964943504052*^9}, {
   3.593965146367446*^9, 3.59396515387148*^9}, {3.593965632781032*^9, 
   3.593965658412372*^9}, {3.593965717252866*^9, 3.593965765732027*^9}, {
   3.593965989042667*^9, 3.59396603681874*^9}, 3.593966944581818*^9, {
   3.594040384945611*^9, 3.594040582644012*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \[NoBreak]\\!\\(380247\\)\[NoBreak] function \
evaluations on \[NoBreak]\\!\\(380\\)\[NoBreak] subregions.\"\>"}]], "Message",\
 "MSG",
 CellChangeTimes->{{3.593966031795422*^9, 3.593966037401582*^9}, 
   3.5939669832105207`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"{", 
   RowBox[{
   "2.6902393005756685`*^-8", ",", "2.686644221390067`*^-11", ",", "1.`"}], 
   "}"}], "}"}]], "Output",
 CellChangeTimes->{
  3.593964952513726*^9, 3.593965012102688*^9, 3.593965147889002*^9, 
   3.5939652236222553`*^9, {3.5939656337696533`*^9, 3.5939656590092583`*^9}, {
   3.593965740262784*^9, 3.593965756608182*^9}, 3.593965798378319*^9, {
   3.593965994607967*^9, 3.5939660374034224`*^9}, 3.5939669832119617`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"B4A", "[", 
  RowBox[{"0.01", ",", "0.1", ",", " ", 
   RowBox[{"-", "0.5"}], ",", " ", "3.9", ",", " ", "0.1"}], "]"}]], "Input",
 CellChangeTimes->{{3.5940421163966637`*^9, 3.594042130474712*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \[NoBreak]\\!\\(458213\\)\[NoBreak] function \
evaluations on \[NoBreak]\\!\\(458\\)\[NoBreak] subregions.\"\>"}]], "Message",\
 "MSG",
 CellChangeTimes->{3.5940432677188168`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"{", 
   RowBox[{
   "2.1391106101580735`*^-7", ",", "2.1370466464643485`*^-10", ",", "1.`"}], 
   "}"}], "}"}]], "Output",
 CellChangeTimes->{3.59404218482903*^9, 3.59404313865376*^9, 
  3.5940432677218723`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"B4E", "[", 
  RowBox[{"0.01", ",", "0.1", ",", " ", 
   RowBox[{"-", "0.5"}], ",", " ", "3.9", ",", " ", "0.1", ",", " ", "0.9", 
   ",", " ", "3.1"}], "]"}]], "Input",
 CellChangeTimes->{{3.5939650317637*^9, 3.5939650658798018`*^9}, {
   3.593965193231915*^9, 3.5939651935674686`*^9}, {3.593965358262787*^9, 
   3.593965367310998*^9}, {3.5939654356857777`*^9, 3.593965437637498*^9}, {
   3.593965711085924*^9, 3.593965711892088*^9}, {3.593965807723896*^9, 
   3.593965838803783*^9}, 3.593965923627893*^9, {3.5939669911341963`*^9, 
   3.593966992813484*^9}, 3.593967059805388*^9, 3.593967138876779*^9, 
   3.594040599939519*^9, {3.594041793251381*^9, 3.594041793846858*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \[NoBreak]\\!\\(100121\\)\[NoBreak] function \
evaluations on \[NoBreak]\\!\\(100\\)\[NoBreak] subregions.\"\>"}]], "Message",\
 "MSG",
 CellChangeTimes->{
  3.593965132508169*^9, 3.593965295003614*^9, 3.5939653600954123`*^9, 
   3.593965425654785*^9, 3.593965536899605*^9, 3.593965712643793*^9, {
   3.5939658156953917`*^9, 3.593965831070055*^9}, 3.5939658926918583`*^9, 
   3.593965977628027*^9, 3.593967053219335*^9, 3.593967127816297*^9, 
   3.5939672091347313`*^9, 3.594041769073805*^9, 3.594041812324151*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "0.18699898472433002`"}], ",", "1.8869665365471202`*^-7"}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.593965132509502*^9, 3.593965295005362*^9, 3.593965360096526*^9, 
   3.593965425656178*^9, 3.593965536900738*^9, 3.593965712645433*^9, {
   3.593965815697122*^9, 3.5939658310716352`*^9}, 3.59396589269307*^9, 
   3.593965977629285*^9, 3.593967053220676*^9, 3.593967127817781*^9, 
   3.5939672091360273`*^9, 3.5940417690765867`*^9, 3.5940418123265257`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"B4ESample1", " ", "=", " ", 
   RowBox[{"ParallelTable", "[", 
    RowBox[{
     RowBox[{"B4E", "[", 
      RowBox[{"0.01", ",", " ", "0.1", ",", " ", 
       RowBox[{"-", "0.5"}], ",", " ", "3.9", ",", " ", "0.1", ",", " ", "mu",
        ",", " ", "phi"}], "]"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"mu", ",", " ", 
       RowBox[{"-", "1."}], ",", " ", "1.", ",", " ", "0.2"}], "}"}], ",", 
     " ", 
     RowBox[{"{", 
      RowBox[{"phi", ",", " ", "0", ",", " ", 
       RowBox[{"2.", "*", "\[Pi]"}], ",", " ", "0.5"}], "}"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.594042147122528*^9, 3.594042151702716*^9}, {
  3.594042268023718*^9, 3.594042326172698*^9}, {3.594042383542235*^9, 
  3.59404243250795*^9}, {3.5940431086003733`*^9, 3.594043123922537*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(448303\\) function evaluations on \
\\!\\(448\\) subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.594043492099641*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(438245\\) function evaluations on \
\\!\\(438\\) subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.594043505489801*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(443227\\) function evaluations on \
\\!\\(443\\) subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.594043512727605*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(485242\\) function evaluations on \
\\!\\(485\\) subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.594043527014786*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(448303\\) function evaluations on \
\\!\\(448\\) subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.594043704250989*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(447219\\) function evaluations on \
\\!\\(447\\) subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.594043767882406*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(541332\\) function evaluations on \
\\!\\(541\\) subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.5940437840191383`*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(501409\\) function evaluations on \
\\!\\(501\\) subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.5940437856081953`*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(448303\\) function evaluations on \
\\!\\(448\\) subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.59404391940825*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"General", "::", "stop"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Further output of \\!\\(\\*StyleBox[\\(Suave :: \
success\\), \\\"MessageName\\\"]\\) will be suppressed during this \
calculation.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.594043919599845*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(471273\\) function evaluations on \
\\!\\(471\\) subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.594044024260138*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"General", "::", "stop"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Further output of \\!\\(\\*StyleBox[\\(Suave :: \
success\\), \\\"MessageName\\\"]\\) will be suppressed during this \
calculation.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.594044024367167*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(484336\\) function evaluations on \
\\!\\(484\\) subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.594044044446166*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"General", "::", "stop"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Further output of \\!\\(\\*StyleBox[\\(Suave :: \
success\\), \\\"MessageName\\\"]\\) will be suppressed during this \
calculation.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.59404404453554*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(494310\\) function evaluations on \
\\!\\(494\\) subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.594044056977088*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"General", "::", "stop"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Further output of \\!\\(\\*StyleBox[\\(Suave :: \
success\\), \\\"MessageName\\\"]\\) will be suppressed during this \
calculation.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.5940440572421913`*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(447315\\) function evaluations on \
\\!\\(447\\) subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.5940471848560953`*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(535297\\) function evaluations on \
\\!\\(535\\) subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.594047480100594*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(573373\\) function evaluations on \
\\!\\(573\\) subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.594047790686908*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"General", "::", "stop"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Further output of \\!\\(\\*StyleBox[\\(Suave :: \
success\\), \\\"MessageName\\\"]\\) will be suppressed during this \
calculation.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.594047790821013*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(590279\\) function evaluations on \
\\!\\(590\\) subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.594050083678689*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(625376\\) function evaluations on \
\\!\\(625\\) subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.5940504356417027`*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(505331\\) function evaluations on \
\\!\\(505\\) subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.594050602172625*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(529300\\) function evaluations on \
\\!\\(529\\) subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.5940506248750887`*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(558252\\) function evaluations on \
\\!\\(558\\) subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.594050735014697*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"General", "::", "stop"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Further output of \\!\\(\\*StyleBox[\\(Suave :: \
success\\), \\\"MessageName\\\"]\\) will be suppressed during this \
calculation.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.594050735228568*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(505331\\) function evaluations on \
\\!\\(505\\) subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.594050822698327*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(485295\\) function evaluations on \
\\!\\(485\\) subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.594050855903891*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(505331\\) function evaluations on \
\\!\\(505\\) subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.594051023183497*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"General", "::", "stop"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Further output of \\!\\(\\*StyleBox[\\(Suave :: \
success\\), \\\"MessageName\\\"]\\) will be suppressed during this \
calculation.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.594051023320961*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(475282\\) function evaluations on \
\\!\\(475\\) subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.594051051361396*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"General", "::", "stop"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Further output of \\!\\(\\*StyleBox[\\(Suave :: \
success\\), \\\"MessageName\\\"]\\) will be suppressed during this \
calculation.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.59405105144587*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"file", " ", "=", " ", 
   RowBox[{
   "OpenWrite", "[", "\"\<~/Research/feedersim/Theory/B4E1.exact\>\"", 
    "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.594042905645515*^9, 3.594042914638812*^9}, {
  3.5940554442303267`*^9, 3.59405544505276*^9}, {3.594074676188671*^9, 
  3.594074724768681*^9}, {3.594075609232737*^9, 3.594075613489286*^9}, {
  3.594075666184463*^9, 3.594075704832814*^9}}],

Cell[BoxData[
 RowBox[{"Write", "[", 
  RowBox[{"file", ",", " ", "B4E1"}], "]"}]], "Input",
 CellChangeTimes->{{3.594075708321197*^9, 3.594075713773102*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"B4E1", " ", "=", " ", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"B4ESample1", "[", 
       RowBox[{"[", "a", "]"}], "]"}], "[", 
      RowBox[{"[", "j", "]"}], "]"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"j", ",", " ", "1", ",", " ", 
       RowBox[{"Length", "[", 
        RowBox[{"B4ESample1", "[", 
         RowBox[{"[", "1", "]"}], "]"}], "]"}]}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"a", ",", " ", "1", ",", " ", 
       RowBox[{"Length", "[", "B4ESample1", "]"}]}], "}"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.594042674604986*^9, 3.594042797056509*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{
   RowBox[{"ListPlot", "[", 
    RowBox[{"B4E1", ",", " ", 
     RowBox[{"AxesLabel", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
       "\"\<Cos(\!\(\*SubscriptBox[\(\[Theta]\), \(12\)]\))\>\"", ",", " ", 
        "\"\<Integral\>\""}], "}"}]}], ",", " ", 
     RowBox[{
     "PlotLabel", "\[Rule]", 
      "\"\<Approx vs Exact (D.9), (k=0.01, p1=0.1, p1\[Mu]=-0.5, \
p1\[Phi]=3.9, p2=0.1\>\""}]}], "]"}], ",", " ", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"2.13911", " ", 
      SuperscriptBox["10", 
       RowBox[{"-", "7"}]]}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"x", " ", ",", " ", 
       RowBox[{"-", "1"}], ",", " ", "1"}], "}"}]}], "]"}]}], "]"}]], "Input",\

 CellChangeTimes->{{3.594042798469531*^9, 3.594042814614847*^9}, {
  3.594042853665738*^9, 3.594042920643702*^9}, {3.594063146455011*^9, 
  3.594063176181329*^9}, {3.594074070535865*^9, 3.594074187605446*^9}, {
  3.59407486617006*^9, 3.59407491018056*^9}, {3.594074961287222*^9, 
  3.594075017089809*^9}}],

Cell[BoxData[
 GraphicsBox[{{{}, {
     {RGBColor[0.24720000000000014`, 0.24, 0.6], 
      PointBox[{{0.5, 2.4953599090080287`*^-7}, {0.022794509880943448`, 
       2.424970312593162*^-7}, {-0.20294065349207552`, 
       2.3563461068746207`*^-7}, {-0.37619090367449814`, 
       2.288769511522829*^-7}, {-0.5159739859787538, 
       2.2256996517887045`*^-7}, {-0.6286758168650943, 
       2.1676617516906452`*^-7}, {-0.7159739859787538, 
       2.114037575738299*^-7}, {-0.7761909036744981, 
       2.0598157309965383`*^-7}, {-0.8029406534920754, 
       2.0112529355203985`*^-7}, {-0.7772054901190566, 
       1.953515570287556*^-7}, {-0.5, 1.8858588550943394`*^-7}}]}, 
     {RGBColor[0.6, 0.24, 0.4428931686004542], 
      PointBox[{{0.5, 2.4953599090080287`*^-7}, {-0.10236307706401593`, 
       2.4329196502342526`*^-7}, {-0.3698174360853546, 
       2.3663264763409566`*^-7}, {-0.5673722756655959, 
       2.3008634491339776`*^-7}, {-0.7203554696142016, 
       2.243255629703394*^-7}, {-0.8372717951066933, 
       2.1977323847354174`*^-7}, {-0.9203554696142017, 
       2.1804987282972217`*^-7}, {-0.9673722756655957, 
       2.2690389973011067`*^-7}, {-0.9698174360853545, 
       2.2276715977371655`*^-7}, {-0.902363077064016, 
       1.9949125374447866`*^-7}, {-0.5, 1.8858588550943394`*^-7}}]}, 
     {RGBColor[0.6, 0.5470136627990908, 0.24], 
      PointBox[{{0.5, 2.4953599090080287`*^-7}, {-0.10452466221888301`, 
       2.433701154171308*^-7}, {-0.37269954962517743`, 
       2.3633655157980257`*^-7}, {-0.5706741515300197, 
       2.2996701457205276`*^-7}, {-0.7238853233908679, 
       2.2428680647879184`*^-7}, {-0.8408744370314718, 
       2.1981826344890037`*^-7}, {-0.9238853233908679, 
       2.1835844464654094`*^-7}, {-0.9706741515300195, 
       2.2980944299248855`*^-7}, {-0.9726995496251774, 
       2.2593372059173156`*^-7}, {-0.9045246622188831, 
       1.9952500683652109`*^-7}, {-0.5, 1.8858588550943394`*^-7}}]}, 
     {RGBColor[0.24, 0.6, 0.33692049419863584`], 
      PointBox[{{0.5, 2.4953599090080287`*^-7}, {0.016838985850171295`, 
       2.426944627527236*^-7}, {-0.21088135219977172`, 
       2.3585936218129336`*^-7}, {-0.38528811689921216`, 
       2.2907503781873583`*^-7}, {-0.5256993159962705, 
       2.2262024580022525`*^-7}, {-0.6386016902497146, 
       2.1689714799743216`*^-7}, {-0.7256993159962706, 
       2.1138187525989645`*^-7}, {-0.785288116899212, 
       2.0629811210176653`*^-7}, {-0.8108813521997716, 
       2.008141825932109*^-7}, {-0.7831610141498287, 
       1.9560540076829573`*^-7}, {-0.5, 1.8858588550943394`*^-7}}]}, 
     {RGBColor[0.24, 0.3531726744018182, 0.6], 
      PointBox[{{0.5, 2.4953599090080287`*^-7}, {0.232013813390643, 
       2.418546929816831*^-7}, {0.07601841785419061, 
       2.3461977802710623`*^-7}, {-0.0566031386147853, 
       2.2765639254530626`*^-7}, {-0.17432029401925397`, 
       2.2108869871960988`*^-7}, {-0.2799769776822617, 
       2.1503867352356507`*^-7}, {-0.37432029401925404`, 
       2.0913483318222797`*^-7}, {-0.4566031386147853, 
       2.0385763763274156`*^-7}, {-0.5239815821458094, 
       1.9846784486881444`*^-7}, {-0.5679861866093571, 
       1.939092648238561*^-7}, {-0.5, 1.8858588550943394`*^-7}}]}, 
     {RGBColor[0.6, 0.24, 0.5632658430022722], 
      PointBox[{{0.5, 2.4953599090080287`*^-7}, {0.48831751813616114`, 
       2.4097546328341536`*^-7}, {0.4177566908482148, 
       2.333274264030903*^-7}, {0.33490723734986455`, 
       2.260895312790918*^-7}, {0.24422190318839598`, 
       2.1951799036028603`*^-7}, {0.1471958635602685, 
       2.1329756234229833`*^-7}, {0.04422190318839592, 
       2.075117600284798*^-7}, {-0.06509276265013555, 
       2.0244009879492808`*^-7}, {-0.18224330915178522`, 
       1.9718429344293442`*^-7}, {-0.3116824818638389, 
       1.927785395138187*^-7}, {-0.5, 1.8858588550943394`*^-7}}]}, 
     {RGBColor[0.6, 0.42664098839727194`, 0.24], 
      PointBox[{{0.5, 2.4953599090080287`*^-7}, {0.7229980142608206, 
       2.400722871008715*^-7}, {0.7306640190144275, 
       2.3216578916865754`*^-7}, {0.6933876165568789, 
       2.2506564300040969`*^-7}, {0.6274535485808098, 
       2.184459648110698*^-7}, {0.5383300237680343, 
       2.1216690523272815`*^-7}, {0.42745354858080975`, 
       2.0646443380302866`*^-7}, {0.2933876165568786, 
       2.0111198676048512`*^-7}, {0.13066401901442737`, 
       1.964851192896515*^-7}, {-0.07700198573917943, 
       1.9181291711955308`*^-7}, {-0.5, 1.8858588550943394`*^-7}}]}, 
     {RGBColor[0.2634521802031821, 0.6, 0.24], 
      PointBox[{{0.5, 2.4953599090080287`*^-7}, {0.878597331544867, 
       2.396820662882098*^-7}, {0.9381297753931559, 2.315322488760117*^-7}, {
       0.9310694997360766, 2.242612899466501*^-7}, {0.8815461696950344, 
       2.1769385097841083`*^-7}, {0.797662219241445, 2.114609662105346*^-7}, {
       0.6815461696950343, 2.0571536219099206`*^-7}, {0.5310694997360763, 
       2.0051739755201032`*^-7}, {0.3381297753931558, 
       1.9578120495926653`*^-7}, {0.07859733154486692, 
       1.9148565246149177`*^-7}, {-0.5, 1.8858588550943394`*^-7}}]}, 
     {RGBColor[0.24, 0.47354534880363613`, 0.6], 
      PointBox[{{0.5, 2.4953599090080287`*^-7}, {0.9170193304012603, 
       2.3939955497284864`*^-7}, {0.9893591072016805, 
       2.3130700814080345`*^-7}, {0.9897600724397196, 
       2.241811842130545*^-7}, {0.9442890310923427, 
       2.1740069139327898`*^-7}, {0.8616988840021005, 
       2.1142936299209888`*^-7}, {0.7442890310923427, 
       2.056539693227529*^-7}, {0.5897600724397194, 
       2.0037627580475114`*^-7}, {0.3893591072016803, 
       1.9565777192210856`*^-7}, {0.11701933040126022`, 
       1.9133228408420715`*^-7}, {-0.5, 1.8858588550943394`*^-7}}]}, 
     {RGBColor[0.5163614825959097, 0.24, 0.6], 
      PointBox[{{0.5, 2.4953599090080287`*^-7}, {0.8288569654958993, 
       2.3975077909655787`*^-7}, {0.8718092873278656, 
       2.3175115333063356`*^-7}, {0.8550898355646744, 
       2.2427454806055337`*^-7}, {0.800320492068883, 
       2.1793753215050783`*^-7}, {0.714761609159832, 
       2.1159819512324976`*^-7}, {0.6003204920688829, 2.06017728105297*^-7}, {
       0.4550898355646742, 2.0076721596423308`*^-7}, {0.27180928732786547`, 
       1.958644996796243*^-7}, {0.028856965495899167`, 
       1.9146043797819776`*^-7}, {-0.5, 1.8858588550943394`*^-7}}]}, 
     {RGBColor[0.6, 0.3062683139954558, 0.24], 
      PointBox[{{0.5, 2.4953599090080287`*^-7}, {0.6356954585275846, 
       2.4042743582976454`*^-7}, {0.6142606113701128, 2.32734774548327*^-7}, {
       0.5600307598866604, 2.2538293522533912`*^-7}, {0.48488907205593085`, 
       2.188132143497772*^-7}, {0.3928257642126408, 
       2.1273225942151694`*^-7}, {0.2848890720559308, 
       2.0702924507041018`*^-7}, {0.16003075988666027`, 
       2.015335810613829*^-7}, {0.014260611370112641`, 
       1.9667589965372972`*^-7}, {-0.1643045414724155, 
       1.9207270568188576`*^-7}, {-0.5, 1.8858588550943394`*^-7}}]}, 
     {RGBColor[0.3838248546049982, 0.6, 0.24], 
      PointBox[{{0.5, 2.4953599090080287`*^-7}, {0.38482748314522824`, 
       2.412601545910629*^-7}, {0.27976997752697097`, 
       2.3387267585010683`*^-7}, {0.17682359767667163`, 
       2.2661238768950707`*^-7}, {0.07522338372802106, 
       2.2034187769306036`*^-7}, {-0.025287528091286336`, 
       2.1399730862570688`*^-7}, {-0.124776616271979, 
       2.0828187204511612`*^-7}, {-0.22317640232332842`, 
       2.0285035088906825`*^-7}, {-0.3202300224730291, 
       1.9784244678956365`*^-7}, {-0.4151725168547718, 
       1.9309727595245315`*^-7}, {-0.5, 1.8858588550943394`*^-7}}]}, 
     {RGBColor[0.24, 0.5939180232054561, 0.6], 
      PointBox[{{0.5, 2.4953599090080287`*^-7}, {0.13767426904894448`, 
       2.4227767988889974`*^-7}, {-0.04976764126807409, 
       2.3487759682650204`*^-7}, {-0.20070917293927745`, 
       2.281995894820871*^-7}, {-0.32837612482180684`, 
       2.219218909045395*^-7}, {-0.43720955158509256`, 
       2.157220146152022*^-7}, {-0.5283761248218068, 
       2.0969555684182066`*^-7}, {-0.6007091729392774, 
       2.0472210647568776`*^-7}, {-0.6497676412680741, 
       1.9917671956066985`*^-7}, {-0.6623257309510555, 
       1.9382607741856437`*^-7}, {-0.5, 
       1.8858588550943394`*^-7}}]}}, {}}, {{}, {}, 
    {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJxTTMoPSmViYGAwAWIQrWcv9ur///f7XX995H27tccuqP3sauN4BL/8fEtu
+r53cP6VK9yvDia+hfNlD+SLmMu8gfO3Zh9h5b/7Cs5P6DH8cLPpJZz/Y1J1
x377F3D+Yr57kjMePoPz53gZFl3ufArnN+UkVmh6PoHzq5lXlRe9fQTnJzac
V7m24iGcP3Xv0xlzsh7A+bvua/CsUboH5ysLODiFzLsN52/S6sh69PgGnH//
6EYFM7NrcP5TnalG+7ddhvMF61aa6TFehPN1zN7rLP18Bs5vdnl7YWfFcThf
NGTdhEjxQwj7vh3lP6GyC85fevtp78uWdXB+bKDMltvJU+D8nVzVVw+VTraH
8U0fMco1RqyH85U/7P5pdWgnnO+5+OZZv4hDcL6rrxDz163H4fz++OsNBtvP
wPnH/WR5A4QuwvlvFUXrLs26DOeXlmgn9otdg/Ov3xEISjtxA85PqLR5zu5z
G84/rG3edZj1HpyvxrnnxRKPB3C+7J5sV9aOh3A+wy+mlUWPHsH53wP3xh1V
fQLn3zaJl1ZteQrnF4QLCCrdeQbnS055EsEi+wLOXxe9YPnnipdwfrhtvceJ
va/g/I1xuUpWf17D+T+4nJl3hr6F84MuWpvJzngH5ztwaE7g83kP50PzC5wP
AO7ozSo=
      "]]}}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesLabel->{
    FormBox[
    "\"Cos(\\!\\(\\*SubscriptBox[\\(\[Theta]\\), \\(12\\)]\\))\"", 
     TraditionalForm], 
    FormBox["\"Integral\"", TraditionalForm]},
  AxesOrigin->{0., 1.8858588550943394`*^-7},
  ImageSize->{517.3333333333339, Automatic},
  Method->{},
  PlotLabel->FormBox[
   "\"Approx vs Exact (D.9), (k=0.01, p1=0.1, p1\[Mu]=-0.5, p1\[Phi]=3.9, \
p2=0.1\"", TraditionalForm],
  PlotRange->{{-0.9726995496251774, 0.9897600724397196}, {
   1.8858588550943394`*^-7, 2.4953599090080287`*^-7}},
  PlotRangeClipping->True,
  PlotRangePadding->{{0.03924919244129794, 0.03924919244129794}, {
   1.2190021078273788`*^-9, 1.2190021078273788`*^-9}}]], "Output",
 CellChangeTimes->{3.594074188110861*^9, 3.5940749122416697`*^9, 
  3.594075017812325*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Start the import of library", "Subsection",
 CellChangeTimes->{{3.593446023119515*^9, 3.593446038368087*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Install", "[", "\"\<~/Research/feedersim/Theory/Suave\>\"", "]"}],
    "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"ParallelEvaluate", "[", 
     RowBox[{
     "Install", "[", "\"\<~/Research/feedersim/Theory/Suave\>\"", "]"}], 
     "]"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"Install", "[", "\"\<Vegas\>\"", "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"Install", "[", "\"\<Divonne\>\"", "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"Install", "[", "\"\<Cuhre\>\"", "]"}], ";"}], 
    "*)"}]}]}]], "Input",
 CellChangeTimes->{{3.5927596892664843`*^9, 3.5927597078205957`*^9}, {
  3.592769356182355*^9, 3.592769381989523*^9}, {3.593445566992251*^9, 
  3.5934456196034946`*^9}, {3.593448309028468*^9, 3.593448327143906*^9}, {
  3.5934696444623213`*^9, 3.593469659743586*^9}, {3.593533761021052*^9, 
  3.5935337677429543`*^9}, {3.594040405964889*^9, 3.5940404064157457`*^9}, {
  3.594040917442429*^9, 3.594040971883329*^9}, {3.594041191623125*^9, 
  3.594041257143462*^9}, {3.594042499624783*^9, 3.594042504973761*^9}, {
  3.594043189256823*^9, 3.5940431978445253`*^9}, {3.594491390701913*^9, 
  3.594491410064373*^9}, {3.594501991532092*^9, 3.594502000186346*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"R", "=", "8."}], ";", " ", 
  RowBox[{"NRS", "=", "50"}], ";", " ", 
  RowBox[{"NN", "=", "100"}], ";"}]], "Input",
 CellChangeTimes->{{3.592771085069105*^9, 3.5927710853591337`*^9}, {
   3.592771129926832*^9, 3.592771154746641*^9}, {3.592771738804964*^9, 
   3.592771740716847*^9}, {3.592773562041222*^9, 3.5927735624464083`*^9}, {
   3.5928275067837563`*^9, 3.592827509507473*^9}, 3.5934443537436857`*^9, {
   3.593448531119543*^9, 3.5934485314273376`*^9}, {3.593469476704803*^9, 
   3.593469477663097*^9}, {3.593516432034074*^9, 3.593516432297584*^9}, {
   3.594501868375342*^9, 3.594501878429996*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"s8", "=", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"Suave", "[", 
     RowBox[{
      RowBox[{"ISigSq", "[", 
       RowBox[{"k", ",", "8."}], "]"}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"k", ",", " ", "kmin", ",", " ", "kmax"}], "}"}], ",", " ", 
      RowBox[{"Verbose", "\[Rule]", "0"}]}], "]"}], ",", 
    RowBox[{"Vegas", "[", 
     RowBox[{
      RowBox[{"ISigSq", "[", 
       RowBox[{"k", ",", "8."}], "]"}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"k", ",", " ", "kmin", ",", " ", "kmax"}], "}"}], ",", " ", 
      RowBox[{"Verbose", "\[Rule]", "0"}]}], "]"}], ",", " ", 
    RowBox[{"Divonne", "[", 
     RowBox[{
      RowBox[{"ISigSq", "[", 
       RowBox[{"k", ",", "8."}], "]"}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"k", ",", " ", "kmin", ",", " ", "kmax"}], "}"}], ",", " ", 
      RowBox[{"Verbose", "\[Rule]", "0"}]}], "]"}], ",", " ", 
    RowBox[{"Cuhre", "[", 
     RowBox[{
      RowBox[{"ISigSq", "[", 
       RowBox[{"k", ",", "8."}], "]"}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"k", ",", " ", "kmin", ",", " ", "kmax"}], "}"}], ",", " ", 
      RowBox[{"Verbose", "\[Rule]", "0"}]}], "]"}]}], "}"}]}]], "Input",
 CellChangeTimes->{{3.5928216804437513`*^9, 3.592821715570931*^9}, {
  3.5928217541675053`*^9, 3.5928218545935717`*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \[NoBreak]\\!\\(10000\\)\[NoBreak] function \
evaluations on \[NoBreak]\\!\\(10\\)\[NoBreak] subregions.\"\>"}]], "Message",\
 "MSG",
 CellChangeTimes->{
  3.592821716967567*^9, {3.592821768939665*^9, 3.592821814181319*^9}, {
   3.592821856031817*^9, 3.592821883629813*^9}, 3.592822038952778*^9, 
   3.592822185363241*^9, 3.592827270492036*^9, 3.593385310152178*^9, 
   3.593387794428018*^9, 3.593446716524247*^9, 3.593469535781733*^9, 
   3.5934696632303133`*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Vegas", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \[NoBreak]\\!\\(27000\\)\[NoBreak] function \
evaluations.\"\>"}]], "Message", "MSG",
 CellChangeTimes->{
  3.592821716967567*^9, {3.592821768939665*^9, 3.592821814181319*^9}, {
   3.592821856031817*^9, 3.592821883629813*^9}, 3.592822038952778*^9, 
   3.592822185363241*^9, 3.592827270492036*^9, 3.593385310152178*^9, 
   3.593387794428018*^9, 3.593446716524247*^9, 3.593469535781733*^9, 
   3.593469665460923*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Divonne", "::", "baddim"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Cannot integrate in \[NoBreak]\\!\\(1\\)\[NoBreak] \
dimensions.\"\>"}]], "Message", "MSG",
 CellChangeTimes->{
  3.592821716967567*^9, {3.592821768939665*^9, 3.592821814181319*^9}, {
   3.592821856031817*^9, 3.592821883629813*^9}, 3.592822038952778*^9, 
   3.592822185363241*^9, 3.592827270492036*^9, 3.593385310152178*^9, 
   3.593387794428018*^9, 3.593446716524247*^9, 3.593469535781733*^9, 
   3.5934696655678368`*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Cuhre", "::", "baddim"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Cannot integrate in \[NoBreak]\\!\\(1\\)\[NoBreak] \
dimensions.\"\>"}]], "Message", "MSG",
 CellChangeTimes->{
  3.592821716967567*^9, {3.592821768939665*^9, 3.592821814181319*^9}, {
   3.592821856031817*^9, 3.592821883629813*^9}, 3.592822038952778*^9, 
   3.592822185363241*^9, 3.592827270492036*^9, 3.593385310152178*^9, 
   3.593387794428018*^9, 3.593446716524247*^9, 3.593469535781733*^9, 
   3.593469665672423*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"{", 
     RowBox[{
     "7.776228260730613`*^8", ",", "736859.5046232337`", ",", 
      "0.5771531696656633`"}], "}"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"{", 
     RowBox[{
     "7.777740004531363`*^8", ",", "723947.9463480129`", ",", 
      "0.0002657882108363705`"}], "}"}], "}"}], ",", "$Failed", ",", 
   "$Failed"}], "}"}]], "Output",
 CellChangeTimes->{{3.592821699845924*^9, 3.592821716970093*^9}, {
   3.592821771854151*^9, 3.592821816966223*^9}, {3.592821859171577*^9, 
   3.592821886763913*^9}, 3.592822042029595*^9, 3.592822188471784*^9, 
   3.59282727451416*^9, 3.5933853126167593`*^9, 3.5933877968528748`*^9, 
   3.59344671920201*^9, 3.5934695360409822`*^9, 3.593469665673601*^9}]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"MP1", "=", "500000"}], ";", 
  RowBox[{"MP2", "=", "10000000"}], ";", " ", 
  RowBox[{"MP3", "=", "100000000"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"IntegratedValue31", "[", 
    RowBox[{"k_", ",", "R_"}], "]"}], ":=", 
   RowBox[{"Suave", "[", 
    RowBox[{
     RowBox[{"I31Value", "[", 
      RowBox[{"k", ",", "R", ",", "q1", ",", " ", "\[Mu]1"}], "]"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"q1", ",", " ", "kmin", ",", " ", "kmax"}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"\[Mu]1", ",", " ", 
       RowBox[{"-", "1."}], ",", "1."}], "}"}], ",", " ", 
     RowBox[{"Verbose", "\[Rule]", "0"}], ",", " ", 
     RowBox[{"MaxPoints", "\[Rule]", "MP1"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"IntegratedValue31D", "[", 
    RowBox[{"k_", ",", "R_"}], "]"}], ":=", 
   RowBox[{"Suave", "[", 
    RowBox[{
     RowBox[{"I31DLog", "[", 
      RowBox[{"k", ",", "R", ",", "q1", ",", " ", "\[Mu]1"}], "]"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"q1", ",", " ", "kmin", ",", " ", "kmax"}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"\[Mu]1", ",", " ", 
       RowBox[{"-", "1."}], ",", "1."}], "}"}], ",", " ", 
     RowBox[{"Verbose", "\[Rule]", "0"}], ",", " ", 
     RowBox[{"MaxPoints", "\[Rule]", "MP1"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"IntegratedValue32", "[", 
    RowBox[{"k_", ",", "R_"}], "]"}], ":=", "\[IndentingNewLine]", 
   RowBox[{"Suave", "[", 
    RowBox[{
     RowBox[{"I32Value", "[", 
      RowBox[{
      "k", ",", "R", ",", "q2", ",", " ", "\[Mu]2", ",", " ", "\[Phi]2", ",", 
       " ", "q1", ",", " ", "\[Mu]1", ",", " ", "\[Phi]1"}], "]"}], ",", " ", 
     
     RowBox[{"{", 
      RowBox[{"q2", ",", " ", "kmin", ",", " ", "kmax"}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"\[Mu]2", ",", " ", 
       RowBox[{"-", "1."}], ",", " ", "1."}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"\[Phi]2", ",", " ", "0.", ",", " ", 
       RowBox[{"2", "*", "\[Pi]"}]}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"q1", ",", " ", "kmin", ",", " ", "kmax"}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"\[Mu]1", ",", " ", 
       RowBox[{"-", "1."}], ",", " ", "1."}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"\[Phi]1", ",", " ", "0.", ",", " ", 
       RowBox[{"2", "*", "\[Pi]"}]}], "}"}], ",", " ", 
     RowBox[{"Verbose", "\[Rule]", "0"}], ",", " ", 
     RowBox[{"MaxPoints", "\[Rule]", "MP2"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"IntegratedValue32D", "[", 
    RowBox[{"k_", ",", "R_"}], "]"}], ":=", "\[IndentingNewLine]", 
   RowBox[{"Suave", "[", 
    RowBox[{
     RowBox[{"I32DLog", "[", 
      RowBox[{
      "k", ",", "R", ",", "q2", ",", " ", "\[Mu]2", ",", " ", "\[Phi]2", ",", 
       " ", "q1", ",", " ", "\[Mu]1", ",", " ", "\[Phi]1"}], "]"}], ",", " ", 
     
     RowBox[{"{", 
      RowBox[{"q2", ",", " ", "kmin", ",", " ", "kmax"}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"\[Mu]2", ",", " ", 
       RowBox[{"-", "1."}], ",", " ", "1."}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"\[Phi]2", ",", " ", "0.", ",", " ", 
       RowBox[{"2", "*", "\[Pi]"}]}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"q1", ",", " ", "kmin", ",", " ", "kmax"}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"\[Mu]1", ",", " ", 
       RowBox[{"-", "1."}], ",", " ", "1."}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"\[Phi]1", ",", " ", "0.", ",", " ", 
       RowBox[{"2", "*", "\[Pi]"}]}], "}"}], ",", " ", 
     RowBox[{"Verbose", "\[Rule]", "0"}], ",", " ", 
     RowBox[{"MaxPoints", "\[Rule]", "MP2"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"IntegratedValue42", "[", 
    RowBox[{"k_", ",", "R_"}], "]"}], ":=", 
   RowBox[{"Suave", "[", 
    RowBox[{
     RowBox[{"I42Value", "[", 
      RowBox[{
      "k", ",", " ", "R", ",", "q2", ",", "\[Mu]2", ",", " ", "\[Phi]2", ",", 
       " ", "q1", ",", " ", "\[Mu]1", ",", " ", "\[Phi]1", ",", " ", "p", ",",
        " ", "\[Mu]", ",", " ", "\[Phi]"}], "]"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"q2", ",", " ", "kmin", ",", " ", "kmax"}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"\[Mu]2", ",", " ", 
       RowBox[{"-", "1."}], ",", "1."}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"\[Phi]2", ",", " ", "0.", ",", " ", 
       RowBox[{"2", "*", "\[Pi]"}]}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"q1", ",", " ", "kmin", ",", " ", "kmax"}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"\[Mu]1", ",", " ", 
       RowBox[{"-", "1."}], ",", " ", "1."}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"\[Phi]1", ",", " ", "0.", ",", " ", 
       RowBox[{"2", "*", "\[Pi]"}]}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"p", ",", " ", "kmin", ",", " ", "kmax"}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"\[Mu]", ",", " ", 
       RowBox[{"-", "1."}], ",", " ", "1."}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"\[Phi]", ",", " ", "0.", ",", " ", 
       RowBox[{"2", "*", "\[Pi]"}]}], "}"}], ",", " ", 
     RowBox[{"Verbose", "\[Rule]", "0"}], ",", " ", 
     RowBox[{"MaxPoints", "\[Rule]", "1000000000"}]}], "]"}]}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.5927713822587433`*^9, 3.592771428110345*^9}, {
   3.59277178643764*^9, 3.5927718002163486`*^9}, 3.592771835166047*^9, {
   3.592771922735601*^9, 3.592771965231061*^9}, 3.592772005506607*^9, {
   3.5927720564287577`*^9, 3.5927720575354567`*^9}, {3.5927721313042927`*^9, 
   3.5927721448949213`*^9}, 3.592772246139324*^9, 3.592772291811881*^9, 
   3.59277276299395*^9, {3.5927727986842012`*^9, 3.5927728156102123`*^9}, 
   3.5927729854591637`*^9, {3.5927732013324203`*^9, 3.59277320148761*^9}, {
   3.5927764813170223`*^9, 3.5927765426969023`*^9}, 3.592776636951676*^9, {
   3.5927804716152287`*^9, 3.592780472541325*^9}, {3.5927927403913507`*^9, 
   3.592792741350944*^9}, {3.5928274890648193`*^9, 3.5928274894709673`*^9}, {
   3.5928698194572563`*^9, 3.592869856128882*^9}, {3.5928716928314457`*^9, 
   3.592871709729251*^9}, 3.5929214998189087`*^9, {3.5929219021096888`*^9, 
   3.592921902492672*^9}, {3.592925999683341*^9, 3.592926000292861*^9}, {
   3.5929265564413567`*^9, 3.592926556831102*^9}, {3.593191581055113*^9, 
   3.593191581475861*^9}, {3.593191754533781*^9, 3.5931917830645647`*^9}, {
   3.593191815502367*^9, 3.593191816072352*^9}, {3.593191878655615*^9, 
   3.593191878791275*^9}, {3.5931924976555853`*^9, 3.593192498517729*^9}, {
   3.593192586056497*^9, 3.5931926404420357`*^9}, {3.593192671597899*^9, 
   3.593192681269909*^9}, {3.593192711929377*^9, 3.593192757830921*^9}, {
   3.593198659354298*^9, 3.59319865969238*^9}, {3.593198760192009*^9, 
   3.593198773002982*^9}, {3.5932145538572817`*^9, 3.5932145586110497`*^9}, 
   3.593276527361408*^9, 3.593385474068346*^9, {3.593387365605975*^9, 
   3.593387376576971*^9}, {3.593387411137824*^9, 3.5933874213499537`*^9}, {
   3.5933877258246813`*^9, 3.5933877261935596`*^9}, {3.593387776242445*^9, 
   3.593387785582316*^9}, {3.5934440797206707`*^9, 3.5934441304232264`*^9}, 
   3.593444162549739*^9, {3.593445076791622*^9, 3.593445077400723*^9}, {
   3.59344567656181*^9, 3.5934456921088467`*^9}, {3.593445743740451*^9, 
   3.593445772832899*^9}, {3.593446652995222*^9, 3.5934466538996277`*^9}, 
   3.593448463917787*^9, {3.593448560945092*^9, 3.59344856217062*^9}, {
   3.593469470777025*^9, 3.5934694716367064`*^9}, {3.593469704675931*^9, 
   3.593469705282894*^9}, 3.593469958821897*^9, {3.593470887233375*^9, 
   3.5934708896625223`*^9}, 3.593472410209291*^9, {3.593515681283361*^9, 
   3.593515690973024*^9}, {3.593515754061695*^9, 3.593515757494254*^9}, {
   3.593515810601774*^9, 3.593515812882571*^9}, {3.593516447686655*^9, 
   3.593516448583859*^9}, {3.593516505426827*^9, 3.593516506605672*^9}, {
   3.593525847869548*^9, 3.593525848351383*^9}, {3.593525916386649*^9, 
   3.59352591844496*^9}, 3.593533741502265*^9, {3.59353396244345*^9, 
   3.593533963949832*^9}, {3.593534002149434*^9, 3.5935340041098013`*^9}, {
   3.593534038350915*^9, 3.593534061441216*^9}, {3.594491943547038*^9, 
   3.594491944617548*^9}, 3.594500956289328*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"12292", "/", "74423"}], " ", "//", "N"}]], "Input",
 CellChangeTimes->{{3.594501302691285*^9, 3.5945013312226887`*^9}, {
  3.5945013957194977`*^9, 3.594501425532916*^9}, {3.594501479044799*^9, 
  3.594501484123478*^9}}],

Cell[BoxData["0.16516399500154522`"], "Output",
 CellChangeTimes->{{3.594501311970312*^9, 3.594501331659449*^9}, {
  3.59450139931439*^9, 3.5945014260634003`*^9}, {3.5945014818473797`*^9, 
  3.594501484663878*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"IntegratedValue31", "[", 
   RowBox[{"0.01", ",", " ", "8."}], "]"}], " ", "//", "Timing"}]], "Input",
 CellChangeTimes->{{3.5945009489883757`*^9, 3.594500963724926*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \[NoBreak]\\!\\(13000\\)\[NoBreak] function \
evaluations on \[NoBreak]\\!\\(13\\)\[NoBreak] subregions.\"\>"}]], "Message",\
 "MSG",
 CellChangeTimes->{3.594500967616047*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"3.63099999999999978328446559316944330931`6.580626162192517", ",", 
   RowBox[{"{", 
    RowBox[{"{", 
     RowBox[{
     "74423.59687992348`", ",", "68.54672475703866`", ",", 
      "0.9607180217248553`"}], "}"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.594500967618188*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"IntegratedValue31D", "[", 
   RowBox[{"0.01", ",", " ", "8."}], "]"}], " ", "//", "Timing"}]], "Input",
 CellChangeTimes->{{3.594500966273291*^9, 3.5945009760125427`*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \[NoBreak]\\!\\(15000\\)\[NoBreak] function \
evaluations on \[NoBreak]\\!\\(15\\)\[NoBreak] subregions.\"\>"}]], "Message",\
 "MSG",
 CellChangeTimes->{3.594500981291099*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"4.93900000000000005684341886080801486969`6.714238939441172", ",", 
   RowBox[{"{", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", "12291.568412744298`"}], ",", "11.2037384206436`", ",", 
      "1.3726015702709473`*^-9"}], "}"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.594500981293213*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"IntegratedValue32", "[", 
   RowBox[{"0.01", ",", "8."}], "]"}], " ", "//", "Timing"}]], "Input",
 CellChangeTimes->{{3.593515699089589*^9, 3.593515706206738*^9}, {
  3.593533747907609*^9, 3.5935337530225487`*^9}, {3.5945005374374866`*^9, 
  3.5945005450798903`*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \[NoBreak]\\!\\(548213\\)\[NoBreak] function \
evaluations on \[NoBreak]\\!\\(548\\)\[NoBreak] subregions.\"\>"}]], "Message",\
 "MSG",
 CellChangeTimes->{3.593515741217267*^9, 3.593515802464931*^9, 
  3.5935158507008753`*^9, 3.593533833447904*^9, 3.593533984075329*^9, 
  3.5935340287863083`*^9, 3.594500723852669*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"175.19800000000000750333128962665796279907`8.264129057383824", ",", 
   RowBox[{"{", 
    RowBox[{"{", 
     RowBox[{
     "3.8100931539380234`", ",", "0.0038089520746813384`", ",", "1.`"}], 
     "}"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.593515741218677*^9, 3.59351580246635*^9, 
  3.5935158507021723`*^9, 3.5935338334574127`*^9, 3.593533984083725*^9, 
  3.593534028801702*^9, 3.594500723973275*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"IntegratedValue32D", "[", 
   RowBox[{"0.01", ",", " ", "8."}], "]"}], " ", "//", " ", 
  "Timing"}]], "Input",
 CellChangeTimes->{{3.594500555562089*^9, 3.594500584829274*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \[NoBreak]\\!\\(438099\\)\[NoBreak] function \
evaluations on \[NoBreak]\\!\\(438\\)\[NoBreak] subregions.\"\>"}]], "Message",\
 "MSG",
 CellChangeTimes->{3.5945008871463623`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"162.3400000000000034106051316484808921814`8.23102545491071", ",", 
   RowBox[{"{", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", "0.6291151626279162`"}], ",", "0.0006283745877405456`", 
      ",", "1.`"}], "}"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.594500887148703*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"klist", " ", "=", " ", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"N", "[", "kmin", "]"}], "+", 
       RowBox[{
        FractionBox[
         RowBox[{"N", "[", "kmin", "]"}], "2"], "*", 
        RowBox[{"(", 
         RowBox[{"i", "-", "1"}], ")"}]}]}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"i", ",", " ", "1", ",", " ", "NRS"}], "}"}]}], "]"}]}], ";"}],
   " ", "//", "N"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"data1", " ", "=", " ", 
     RowBox[{"ParallelTable", "[", 
      RowBox[{
       RowBox[{"IntegratedValue32", "[", 
        RowBox[{
         RowBox[{"klist", "[", 
          RowBox[{"[", "i", "]"}], "]"}], ",", "R"}], "]"}], ",", " ", 
       RowBox[{"{", 
        RowBox[{"i", ",", " ", "1", ",", " ", "NRS"}], "}"}]}], "]"}]}], 
    ";"}], " ", "//", "AbsoluteTiming"}], "\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{
   "now", " ", "save", " ", "the", " ", "data", " ", "in", " ", "files"}], 
   " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fname", " ", "=", " ", 
   RowBox[{"StringJoin", "[", 
    RowBox[{
     RowBox[{"Directory", "[", "]"}], ",", "\"\</F3s/f32_\>\"", ",", 
     RowBox[{"ToString", "[", "NN", "]"}], ",", "\"\<.txt\>\""}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"file", " ", "=", " ", 
   RowBox[{"OpenWrite", "[", "fname", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Do", "[", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"dataI", "=", 
      RowBox[{"Flatten", "[", 
       RowBox[{"data1", "[", 
        RowBox[{"[", "i", "]"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"WriteString", "[", 
      RowBox[{"file", ",", " ", 
       RowBox[{"StringJoin", "[", 
        RowBox[{
         RowBox[{"ToString", "[", 
          RowBox[{"FortranForm", "[", 
           RowBox[{"klist", "[", 
            RowBox[{"[", "i", "]"}], "]"}], "]"}], "]"}], ",", "\"\< \>\"", 
         ",", 
         RowBox[{"ToString", "[", 
          RowBox[{"FortranForm", "[", 
           RowBox[{"dataI", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "]"}], "]"}], ",", "\"\< \>\"", 
         ",", " ", 
         RowBox[{"ToString", "[", 
          RowBox[{"FortranForm", "[", 
           RowBox[{"dataI", "[", 
            RowBox[{"[", "2", "]"}], "]"}], "]"}], "]"}], ",", " ", 
         "\"\<\\n\>\""}], "]"}]}], "]"}]}], ",", " ", 
    RowBox[{"{", 
     RowBox[{"i", ",", " ", "1", ",", " ", "NRS"}], "}"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Close", "[", "fname", "]"}], 
  ";"}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.5927711753835773`*^9, 3.592771340892243*^9}, {
   3.592771439632832*^9, 3.5927715236642447`*^9}, {3.592771658103341*^9, 
   3.5927717445515823`*^9}, 3.5927719474711227`*^9, {3.592772028051697*^9, 
   3.592772038227812*^9}, {3.592772294916617*^9, 3.592772298437229*^9}, {
   3.592772540372761*^9, 3.5927727604843397`*^9}, {3.5927727922285223`*^9, 
   3.592772794261052*^9}, {3.5927728264783163`*^9, 3.592772840905918*^9}, {
   3.592772961968841*^9, 3.59277300715025*^9}, {3.592773120643577*^9, 
   3.592773168335147*^9}, {3.592773206143713*^9, 3.5927732073818808`*^9}, {
   3.5927735647265253`*^9, 3.592773566210203*^9}, {3.592828381830942*^9, 
   3.592828382160915*^9}, 3.592830825913354*^9, {3.5932145667468987`*^9, 
   3.593214580133106*^9}, {3.593385378589635*^9, 3.593385388859503*^9}, {
   3.593387660286862*^9, 3.593387661045257*^9}, {3.5933876957542877`*^9, 
   3.5933876975308733`*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \[NoBreak]\\!\\(100152\\)\[NoBreak] function \
evaluations on \[NoBreak]\\!\\(100\\)\[NoBreak] subregions.\"\>"}]], "Message",\
 "MSG",
 CellChangeTimes->{
  3.592772134357877*^9, 3.5927721852904797`*^9, 3.592772257972266*^9, 
   3.592772341839302*^9, 3.59277276910349*^9, {3.592772800861311*^9, 
   3.592772820966248*^9}, {3.592772987813488*^9, 3.592773015546142*^9}, {
   3.592773124734326*^9, 3.592773169203025*^9}, 3.59277325710047*^9, 
   3.59277362700327*^9, 3.592828409524033*^9, 3.593214581715234*^9, 
   3.5932146347119513`*^9, {3.593387707578248*^9, 3.59338773180692*^9}, 
   3.5933878236215982`*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \[NoBreak]\\!\\(100179\\)\[NoBreak] function \
evaluations on \[NoBreak]\\!\\(100\\)\[NoBreak] subregions.\"\>"}]], "Message",\
 "MSG",
 CellChangeTimes->{
  3.592772134357877*^9, 3.5927721852904797`*^9, 3.592772257972266*^9, 
   3.592772341839302*^9, 3.59277276910349*^9, {3.592772800861311*^9, 
   3.592772820966248*^9}, {3.592772987813488*^9, 3.592773015546142*^9}, {
   3.592773124734326*^9, 3.592773169203025*^9}, 3.59277325710047*^9, 
   3.59277362700327*^9, 3.592828409524033*^9, 3.593214581715234*^9, 
   3.5932146347119513`*^9, {3.593387707578248*^9, 3.59338773180692*^9}, 
   3.593387850258992*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"klist", " ", "=", " ", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"N", "[", "kmin", "]"}], "+", 
       RowBox[{
        FractionBox[
         RowBox[{"N", "[", "kmin", "]"}], "2"], "*", 
        RowBox[{"(", 
         RowBox[{"i", "-", "1"}], ")"}]}]}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"i", ",", " ", "1", ",", " ", "NRS"}], "}"}]}], "]"}]}], ";"}],
   " ", "//", "N"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"data1", " ", "=", " ", 
     RowBox[{"ParallelTable", "[", 
      RowBox[{
       RowBox[{"IntegratedValue31", "[", 
        RowBox[{
         RowBox[{"klist", "[", 
          RowBox[{"[", "i", "]"}], "]"}], ",", "R"}], "]"}], ",", " ", 
       RowBox[{"{", 
        RowBox[{"i", ",", " ", "1", ",", " ", "NRS"}], "}"}]}], "]"}]}], 
    ";"}], " ", "//", "AbsoluteTiming"}], "\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{
   "now", " ", "save", " ", "the", " ", "data", " ", "in", " ", "files"}], 
   " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fname", " ", "=", " ", 
   RowBox[{"StringJoin", "[", 
    RowBox[{
     RowBox[{"Directory", "[", "]"}], ",", "\"\</F3s/f31_\>\"", ",", 
     RowBox[{"ToString", "[", "NN", "]"}], ",", "\"\<.txt\>\""}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"file", " ", "=", " ", 
   RowBox[{"OpenWrite", "[", "fname", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Do", "[", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"dataI", "=", 
      RowBox[{"Flatten", "[", 
       RowBox[{"data1", "[", 
        RowBox[{"[", "i", "]"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"WriteString", "[", 
      RowBox[{"file", ",", " ", 
       RowBox[{"StringJoin", "[", 
        RowBox[{
         RowBox[{"ToString", "[", 
          RowBox[{"FortranForm", "[", 
           RowBox[{"klist", "[", 
            RowBox[{"[", "i", "]"}], "]"}], "]"}], "]"}], ",", "\"\< \>\"", 
         ",", 
         RowBox[{"ToString", "[", 
          RowBox[{"FortranForm", "[", 
           RowBox[{"dataI", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "]"}], "]"}], ",", "\"\< \>\"", 
         ",", " ", 
         RowBox[{"ToString", "[", 
          RowBox[{"FortranForm", "[", 
           RowBox[{"dataI", "[", 
            RowBox[{"[", "2", "]"}], "]"}], "]"}], "]"}], ",", " ", 
         "\"\<\\n\>\""}], "]"}]}], "]"}]}], ",", " ", 
    RowBox[{"{", 
     RowBox[{"i", ",", " ", "1", ",", " ", "NRS"}], "}"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Close", "[", "fname", "]"}], ";"}]}], "Input",
 CellChangeTimes->{{3.5927760053853397`*^9, 3.592776007089756*^9}, {
   3.592776555470686*^9, 3.5927765690179167`*^9}, {3.5928275186677313`*^9, 
   3.5928275188262653`*^9}, {3.593385447641788*^9, 3.5933854571511497`*^9}, {
   3.593444178432922*^9, 3.593444181129621*^9}, 3.593444346163642*^9, {
   3.593444378932732*^9, 3.5934444439247427`*^9}, {3.5934451326830263`*^9, 
   3.593445135996913*^9}, {3.5934451933057756`*^9, 3.593445226702444*^9}, {
   3.5934452718116627`*^9, 3.593445274389716*^9}, {3.593445336141431*^9, 
   3.59344533856743*^9}, {3.593445714356337*^9, 3.593445715635779*^9}, {
   3.593445871687708*^9, 3.5934458834966288`*^9}, {3.593446203132511*^9, 
   3.593446203557993*^9}, 3.59344630713479*^9, {3.593446501096034*^9, 
   3.59344657154554*^9}, {3.5934466611365*^9, 3.593446696251484*^9}, 
   3.5934467461033573`*^9, {3.5934467951465883`*^9, 3.59344682621557*^9}, {
   3.593446898006463*^9, 3.593446942161002*^9}, {3.593448387763578*^9, 
   3.593448397601325*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(10000\\) function evaluations on \\!\\(10\\) \
subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.593448344207114*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], ":", 
  " ", "\<\"Needed \\!\\(10000\\) function evaluations on \\!\\(10\\) \
subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.593448344310714*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(10000\\) function evaluations on \\!\\(10\\) \
subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.593448344323024*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(10000\\) function evaluations on \\!\\(10\\) \
subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.5934483445107527`*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(10000\\) function evaluations on \\!\\(10\\) \
subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.5934483445436487`*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(10000\\) function evaluations on \\!\\(10\\) \
subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.593448344887891*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(10000\\) function evaluations on \\!\\(10\\) \
subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.593448346416976*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(10000\\) function evaluations on \\!\\(10\\) \
subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.593448347088873*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(10000\\) function evaluations on \\!\\(10\\) \
subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.59344834782476*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \\!\\(10000\\) function evaluations on \\!\\(10\\) \
subregions.\"\>"}]], "Message", "MSG",
 ShowCellLabel->True,
 CellChangeTimes->{3.593448348578274*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "6.60006700000000012806822269340045750141`6.840148257546129", ",", "Null"}],
   "}"}]], "Output",
 CellChangeTimes->{
  3.592776007295603*^9, 3.592776633380411*^9, 3.592776998575511*^9, 
   3.592793106765875*^9, 3.592827888014913*^9, 3.5933855471433277`*^9, 
   3.593385651621048*^9, 3.593444242919923*^9, {3.5934443737465076`*^9, 
   3.593444435000843*^9}, 3.593444744965992*^9, 3.593445117532407*^9, 
   3.5934451596422987`*^9, 3.593445217738965*^9, 3.593445252600267*^9, 
   3.593445297969771*^9, 3.593445362615986*^9, 3.593445724517857*^9, 
   3.5934457549547873`*^9, 3.593445802151476*^9, 3.593445908032243*^9, 
   3.5934462081657143`*^9, 3.593446271341393*^9, {3.593446513136568*^9, 
   3.5934465746079397`*^9}, {3.593446684471424*^9, 3.59344669987904*^9}, 
   3.593446772753553*^9, 3.593446804919759*^9, 3.5934469201707983`*^9, 
   3.593448349354622*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"klist", " ", "=", " ", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"N", "[", "kmin", "]"}], "+", 
       RowBox[{
        FractionBox[
         RowBox[{"N", "[", "kmin", "]"}], "2"], "*", 
        RowBox[{"(", 
         RowBox[{"i", "-", "1"}], ")"}]}]}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"i", ",", " ", "1", ",", " ", "NRS"}], "}"}]}], "]"}]}], ";"}],
   " ", "//", "N"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"data1", " ", "=", " ", 
     RowBox[{"ParallelTable", "[", 
      RowBox[{
       RowBox[{"IntegratedValue31D", "[", 
        RowBox[{
         RowBox[{"klist", "[", 
          RowBox[{"[", "i", "]"}], "]"}], ",", "R"}], "]"}], ",", " ", 
       RowBox[{"{", 
        RowBox[{"i", ",", " ", "1", ",", " ", "NRS"}], "}"}]}], "]"}]}], 
    ";"}], " ", "//", "AbsoluteTiming"}], "\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{
   "now", " ", "save", " ", "the", " ", "data", " ", "in", " ", "files"}], 
   " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fname", " ", "=", " ", 
   RowBox[{"StringJoin", "[", 
    RowBox[{
     RowBox[{"Directory", "[", "]"}], ",", "\"\</F3s/f31D_\>\"", ",", 
     RowBox[{"ToString", "[", "NN", "]"}], ",", "\"\<.txt\>\""}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"file", " ", "=", " ", 
   RowBox[{"OpenWrite", "[", "fname", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Do", "[", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"dataI", "=", 
      RowBox[{"Flatten", "[", 
       RowBox[{"data1", "[", 
        RowBox[{"[", "i", "]"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"WriteString", "[", 
      RowBox[{"file", ",", " ", 
       RowBox[{"StringJoin", "[", 
        RowBox[{
         RowBox[{"ToString", "[", 
          RowBox[{"FortranForm", "[", 
           RowBox[{"klist", "[", 
            RowBox[{"[", "i", "]"}], "]"}], "]"}], "]"}], ",", "\"\< \>\"", 
         ",", 
         RowBox[{"ToString", "[", 
          RowBox[{"FortranForm", "[", 
           RowBox[{"dataI", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "]"}], "]"}], ",", "\"\< \>\"", 
         ",", " ", 
         RowBox[{"ToString", "[", 
          RowBox[{"FortranForm", "[", 
           RowBox[{"dataI", "[", 
            RowBox[{"[", "2", "]"}], "]"}], "]"}], "]"}], ",", " ", 
         "\"\<\\n\>\""}], "]"}]}], "]"}]}], ",", " ", 
    RowBox[{"{", 
     RowBox[{"i", ",", " ", "1", ",", " ", "NRS"}], "}"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Close", "[", "fname", "]"}], ";"}]}], "Input",
 CellChangeTimes->{{3.592841407925789*^9, 3.592841407951898*^9}, {
  3.5929108243786373`*^9, 3.5929108407572603`*^9}, {3.593387456638483*^9, 
  3.593387464013657*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \[NoBreak]\\!\\(14000\\)\[NoBreak] function \
evaluations on \[NoBreak]\\!\\(14\\)\[NoBreak] subregions.\"\>"}]], "Message",\
 "MSG",
 CellChangeTimes->{3.592914179405549*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \[NoBreak]\\!\\(14000\\)\[NoBreak] function \
evaluations on \[NoBreak]\\!\\(14\\)\[NoBreak] subregions.\"\>"}]], "Message",\
 "MSG",
 CellChangeTimes->{3.5929141843323307`*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Suave", "::", "success"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Needed \[NoBreak]\\!\\(14000\\)\[NoBreak] function \
evaluations on \[NoBreak]\\!\\(14\\)\[NoBreak] subregions.\"\>"}]], "Message",\
 "MSG",
 CellChangeTimes->{3.5929141901592417`*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"General", "::", "stop"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Further output of \[NoBreak]\\!\\(\\*StyleBox[\\(Suave :: \
success\\), \\\"MessageName\\\"]\\)\[NoBreak] will be suppressed during this \
calculation. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \
ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/General/stop\\\", ButtonNote -> \
\\\"General::stop\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.592914191604917*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "412.90184499999998024577507749199867248535`8.636446736748972", ",", 
   "Null"}], "}"}]], "Output",
 CellChangeTimes->{3.592914571553041*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"klist", " ", "=", " ", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"N", "[", "kmin", "]"}], "+", 
       RowBox[{
        FractionBox[
         RowBox[{"N", "[", "kmin", "]"}], "2"], "*", 
        RowBox[{"(", 
         RowBox[{"i", "-", "1"}], ")"}]}]}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"i", ",", " ", "1", ",", " ", "NRS"}], "}"}]}], "]"}]}], ";"}],
   " ", "//", "N"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"data1", " ", "=", " ", 
     RowBox[{"ParallelTable", "[", 
      RowBox[{
       RowBox[{"IntegratedValue32D", "[", 
        RowBox[{
         RowBox[{"klist", "[", 
          RowBox[{"[", "i", "]"}], "]"}], ",", "R"}], "]"}], ",", " ", 
       RowBox[{"{", 
        RowBox[{"i", ",", " ", "1", ",", " ", "NRS"}], "}"}]}], "]"}]}], 
    ";"}], " ", "//", "AbsoluteTiming"}], "\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{
   "now", " ", "save", " ", "the", " ", "data", " ", "in", " ", "files"}], 
   " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fname", " ", "=", " ", 
   RowBox[{"StringJoin", "[", 
    RowBox[{
     RowBox[{"Directory", "[", "]"}], ",", " ", "\"\</F3s/f32D_\>\"", ",", 
     RowBox[{"ToString", "[", "NN", "]"}], ",", "\"\<.txt\>\""}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"file", " ", "=", " ", 
   RowBox[{"OpenWrite", "[", "fname", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Do", "[", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"dataI", "=", 
      RowBox[{"Flatten", "[", 
       RowBox[{"data1", "[", 
        RowBox[{"[", "i", "]"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"WriteString", "[", 
      RowBox[{"file", ",", " ", 
       RowBox[{"StringJoin", "[", 
        RowBox[{
         RowBox[{"ToString", "[", 
          RowBox[{"FortranForm", "[", 
           RowBox[{"klist", "[", 
            RowBox[{"[", "i", "]"}], "]"}], "]"}], "]"}], ",", "\"\< \>\"", 
         ",", 
         RowBox[{"ToString", "[", 
          RowBox[{"FortranForm", "[", 
           RowBox[{"dataI", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "]"}], "]"}], ",", "\"\< \>\"", 
         ",", " ", 
         RowBox[{"ToString", "[", 
          RowBox[{"FortranForm", "[", 
           RowBox[{"dataI", "[", 
            RowBox[{"[", "2", "]"}], "]"}], "]"}], "]"}], ",", " ", 
         "\"\<\\n\>\""}], "]"}]}], "]"}]}], ",", " ", 
    RowBox[{"{", 
     RowBox[{"i", ",", " ", "1", ",", " ", "NRS"}], "}"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Close", "[", "fname", "]"}], ";"}]}], "Input",
 CellChangeTimes->{{3.592910853414578*^9, 3.592910857623952*^9}, {
  3.5933874711779213`*^9, 3.593387480007723*^9}}],

Cell[BoxData["$Aborted"], "Output",
 CellChangeTimes->{3.592918123241426*^9, 3.593222500253215*^9, 
  3.593387486638891*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"klist", " ", "=", " ", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"N", "[", "kmin", "]"}], "+", 
       RowBox[{
        RowBox[{"N", "[", "kmin", "]"}], "*", 
        RowBox[{"(", 
         RowBox[{"i", "-", "1"}], ")"}]}]}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"i", ",", " ", "1", ",", " ", "NRS"}], "}"}]}], "]"}]}], ";"}],
   " ", "//", "N"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"data1", " ", "=", " ", 
     RowBox[{"ParallelTable", "[", 
      RowBox[{
       RowBox[{"IntegratedValue42", "[", 
        RowBox[{
         RowBox[{"klist", "[", 
          RowBox[{"[", "i", "]"}], "]"}], ",", "R"}], "]"}], ",", " ", 
       RowBox[{"{", 
        RowBox[{"i", ",", " ", "1", ",", " ", "NRS"}], "}"}]}], "]"}]}], 
    ";"}], " ", "//", "AbsoluteTiming"}], "\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{
   "now", " ", "save", " ", "the", " ", "data", " ", "in", " ", "files"}], 
   " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fname", " ", "=", " ", 
   RowBox[{"StringJoin", "[", 
    RowBox[{
     RowBox[{"Directory", "[", "]"}], ",", "\"\</F3s/f42_\>\"", ",", 
     RowBox[{"ToString", "[", "NN", "]"}], ",", "\"\<_8.txt\>\""}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"file", " ", "=", " ", 
   RowBox[{"OpenWrite", "[", "fname", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Do", "[", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"dataI", "=", 
      RowBox[{"Flatten", "[", 
       RowBox[{"data1", "[", 
        RowBox[{"[", "i", "]"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"WriteString", "[", 
      RowBox[{"file", ",", " ", 
       RowBox[{"StringJoin", "[", 
        RowBox[{
         RowBox[{"ToString", "[", 
          RowBox[{"FortranForm", "[", 
           RowBox[{"klist", "[", 
            RowBox[{"[", "i", "]"}], "]"}], "]"}], "]"}], ",", "\"\< \>\"", 
         ",", 
         RowBox[{"ToString", "[", 
          RowBox[{"FortranForm", "[", 
           RowBox[{"dataI", "[", 
            RowBox[{"[", "1", "]"}], "]"}], "]"}], "]"}], ",", "\"\< \>\"", 
         ",", " ", 
         RowBox[{"ToString", "[", 
          RowBox[{"FortranForm", "[", 
           RowBox[{"dataI", "[", 
            RowBox[{"[", "2", "]"}], "]"}], "]"}], "]"}], ",", " ", 
         "\"\<\\n\>\""}], "]"}]}], "]"}]}], ",", " ", 
    RowBox[{"{", 
     RowBox[{"i", ",", " ", "1", ",", " ", "NRS"}], "}"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Close", "[", "fname", "]"}], ";"}]}], "Input",
 CellChangeTimes->{{3.5933876394829607`*^9, 3.5933876479406013`*^9}, {
  3.59338768554923*^9, 3.593387687598267*^9}, {3.593469505910149*^9, 
  3.5934695083318243`*^9}, {3.5935258920314903`*^9, 3.5935258925122957`*^9}}],

Cell[BoxData["$Aborted"], "Output",
 CellChangeTimes->{
  3.593469513919517*^9, 3.5934696940417147`*^9, 3.593469798295259*^9, 
   3.593470733720428*^9, 3.5934719526422253`*^9, 3.593480055911484*^9, 
   3.593516616456401*^9, 3.593516689613307*^9, 3.593519198306422*^9, {
   3.593525906743668*^9, 3.5935259361802053`*^9}, 3.593533730980015*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Part", "::", "partd"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Part specification \[NoBreak]\\!\\(data1 \
\[LeftDoubleBracket] 1 \[RightDoubleBracket]\\)\[NoBreak] is longer than \
depth of object. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\
\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/General/partd\\\", ButtonNote -> \
\\\"Part::partd\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{{3.593525906840941*^9, 3.593525936196621*^9}, 
   3.5935337311584463`*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Part", "::", "partd"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Part specification \[NoBreak]\\!\\(data1 \
\[LeftDoubleBracket] 1 \[RightDoubleBracket]\\)\[NoBreak] is longer than \
depth of object. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\
\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/General/partd\\\", ButtonNote -> \
\\\"Part::partd\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{{3.593525906840941*^9, 3.593525936196621*^9}, 
   3.593533731175675*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Part", "::", "partd"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Part specification \[NoBreak]\\!\\(data1 \
\[LeftDoubleBracket] 2 \[RightDoubleBracket]\\)\[NoBreak] is longer than \
depth of object. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\
\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/General/partd\\\", ButtonNote -> \
\\\"Part::partd\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{{3.593525906840941*^9, 3.593525936196621*^9}, 
   3.593533731189147*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"General", "::", "stop"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Further output of \[NoBreak]\\!\\(\\*StyleBox[\\(Part :: \
partd\\), \\\"MessageName\\\"]\\)\[NoBreak] will be suppressed during this \
calculation. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \
ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/General/stop\\\", ButtonNote -> \
\\\"General::stop\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{{3.593525906840941*^9, 3.593525936196621*^9}, 
   3.593533731444008*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"7633", "*", 
  RowBox[{"10", "/", "2"}]}]], "Input",
 CellChangeTimes->{{3.593516462652853*^9, 3.59351646909225*^9}}],

Cell[BoxData["38165"], "Output",
 CellChangeTimes->{3.593516469568511*^9}]
}, Open  ]]
}, Open  ]]
},
WindowSize->{636, 733},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
Magnification:>0.75 Inherited,
FrontEndVersion->"9.0 for Linux x86 (64-bit) (February 7, 2013)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 539, 7, 43, "Subsubsection"],
Cell[1121, 31, 22173, 654, 855, "Input",
 InitializationCell->True],
Cell[CellGroupData[{
Cell[23319, 689, 197, 4, 18, "Input"],
Cell[23519, 695, 580, 8, 70, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[24136, 708, 494, 12, 21, "Input"],
Cell[24633, 722, 26783, 555, 70, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[51453, 1282, 216, 5, 18, "Input"],
Cell[51672, 1289, 680, 15, 70, "Output"]
}, Open  ]]
}, Closed]],
Cell[CellGroupData[{
Cell[52401, 1310, 203, 2, 44, "Subsection"],
Cell[52607, 1314, 26139, 691, 1554, "Input"],
Cell[78749, 2007, 471, 15, 48, "Input"],
Cell[79223, 2024, 562, 16, 42, "Input"],
Cell[CellGroupData[{
Cell[79810, 2044, 847, 23, 76, "Input"],
Cell[80660, 2069, 286, 8, 18, "Message"],
Cell[80949, 2079, 99, 1, 24, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[81085, 2085, 184, 3, 24, "Input"],
Cell[81272, 2090, 105, 1, 24, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[81414, 2096, 130, 2, 24, "Input"],
Cell[81547, 2100, 330, 8, 18, "Message"],
Cell[81880, 2110, 222, 6, 29, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[82139, 2121, 181, 3, 24, "Input"],
Cell[82323, 2126, 330, 8, 18, "Message"],
Cell[82656, 2136, 282, 7, 24, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[82975, 2148, 219, 4, 24, "Input"],
Cell[83197, 2154, 356, 9, 18, "Message"],
Cell[83556, 2165, 248, 6, 24, "Output"]
}, Open  ]],
Cell[83819, 2174, 2487, 61, 144, "Input"],
Cell[CellGroupData[{
Cell[86331, 2239, 959, 22, 76, "Input"],
Cell[87293, 2263, 358, 9, 18, "Message"],
Cell[87654, 2274, 476, 10, 29, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[88167, 2289, 222, 4, 24, "Input"],
Cell[88392, 2295, 308, 8, 18, "Message"],
Cell[88703, 2305, 251, 7, 29, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[88991, 2317, 698, 11, 24, "Input"],
Cell[89692, 2330, 643, 13, 18, "Message"],
Cell[90338, 2345, 514, 10, 29, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[90889, 2360, 818, 19, 42, "Input"],
Cell[91710, 2381, 286, 8, 18, "Message"],
Cell[91999, 2391, 286, 8, 18, "Message"],
Cell[92288, 2401, 286, 8, 18, "Message"],
Cell[92577, 2411, 286, 8, 18, "Message"],
Cell[92866, 2421, 286, 8, 18, "Message"],
Cell[93155, 2431, 286, 8, 18, "Message"],
Cell[93444, 2441, 288, 8, 18, "Message"],
Cell[93735, 2451, 288, 8, 18, "Message"],
Cell[94026, 2461, 285, 8, 18, "Message"],
Cell[94314, 2471, 343, 9, 18, "Message"],
Cell[94660, 2482, 286, 8, 18, "Message"],
Cell[94949, 2492, 343, 9, 18, "Message"],
Cell[95295, 2503, 286, 8, 18, "Message"],
Cell[95584, 2513, 342, 9, 18, "Message"],
Cell[95929, 2524, 286, 8, 18, "Message"],
Cell[96218, 2534, 345, 9, 18, "Message"],
Cell[96566, 2545, 288, 8, 18, "Message"],
Cell[96857, 2555, 286, 8, 18, "Message"],
Cell[97146, 2565, 286, 8, 18, "Message"],
Cell[97435, 2575, 343, 9, 18, "Message"],
Cell[97781, 2586, 286, 8, 18, "Message"],
Cell[98070, 2596, 288, 8, 18, "Message"],
Cell[98361, 2606, 286, 8, 18, "Message"],
Cell[98650, 2616, 288, 8, 18, "Message"],
Cell[98941, 2626, 286, 8, 18, "Message"],
Cell[99230, 2636, 343, 9, 18, "Message"],
Cell[99576, 2647, 286, 8, 18, "Message"],
Cell[99865, 2657, 286, 8, 18, "Message"],
Cell[100154, 2667, 286, 8, 18, "Message"],
Cell[100443, 2677, 343, 9, 18, "Message"],
Cell[100789, 2688, 286, 8, 18, "Message"],
Cell[101078, 2698, 342, 9, 18, "Message"]
}, Open  ]],
Cell[101435, 2710, 431, 9, 24, "Input"],
Cell[101869, 2721, 158, 3, 24, "Input"],
Cell[102030, 2726, 658, 18, 42, "Input"],
Cell[CellGroupData[{
Cell[102713, 2748, 1046, 27, 81, "Input"],
Cell[103762, 2777, 9844, 175, 250, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[113655, 2958, 115, 1, 31, "Subsection"],
Cell[113773, 2961, 1299, 27, 93, "Input"],
Cell[115075, 2990, 640, 11, 24, "Input"],
Cell[CellGroupData[{
Cell[115740, 3005, 1314, 33, 76, "Input"],
Cell[117057, 3040, 593, 13, 18, "Message"],
Cell[117653, 3055, 543, 12, 18, "Message"],
Cell[118199, 3069, 545, 12, 18, "Message"],
Cell[118747, 3083, 541, 12, 18, "Message"],
Cell[119291, 3097, 772, 18, 49, "Output"]
}, Open  ]],
Cell[120078, 3118, 8288, 174, 348, "Input"],
Cell[CellGroupData[{
Cell[128391, 3296, 254, 5, 24, "Input"],
Cell[128648, 3303, 214, 3, 24, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[128899, 3311, 205, 4, 24, "Input"],
Cell[129107, 3317, 304, 8, 18, "Message"],
Cell[129414, 3327, 322, 8, 24, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[129773, 3340, 206, 4, 24, "Input"],
Cell[129982, 3346, 304, 8, 18, "Message"],
Cell[130289, 3356, 342, 8, 29, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[130668, 3369, 302, 6, 24, "Input"],
Cell[130973, 3377, 448, 10, 18, "Message"],
Cell[131424, 3389, 453, 10, 24, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[131914, 3404, 212, 5, 24, "Input"],
Cell[132129, 3411, 308, 8, 18, "Message"],
Cell[132440, 3421, 327, 8, 24, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[132804, 3434, 3607, 88, 255, "Input"],
Cell[136414, 3524, 739, 15, 18, "Message"],
Cell[137156, 3541, 737, 15, 18, "Message"]
}, Open  ]],
Cell[CellGroupData[{
Cell[137930, 3561, 3648, 88, 238, "Input"],
Cell[141581, 3651, 284, 8, 18, "Message"],
Cell[141868, 3661, 271, 7, 18, "Message"],
Cell[142142, 3670, 284, 8, 18, "Message"],
Cell[142429, 3680, 286, 8, 18, "Message"],
Cell[142718, 3690, 286, 8, 18, "Message"],
Cell[143007, 3700, 284, 8, 18, "Message"],
Cell[143294, 3710, 284, 8, 18, "Message"],
Cell[143581, 3720, 284, 8, 18, "Message"],
Cell[143868, 3730, 283, 8, 18, "Message"],
Cell[144154, 3740, 284, 8, 18, "Message"],
Cell[144441, 3750, 903, 16, 24, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[145381, 3771, 2834, 77, 238, "Input"],
Cell[148218, 3850, 304, 8, 18, "Message"],
Cell[148525, 3860, 306, 8, 18, "Message"],
Cell[148834, 3870, 306, 8, 18, "Message"],
Cell[149143, 3880, 531, 11, 18, "Message"],
Cell[149677, 3893, 184, 5, 24, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[149898, 3903, 2788, 76, 238, "Input"],
Cell[152689, 3981, 124, 2, 24, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[152850, 3988, 2857, 76, 212, "Input"],
Cell[155710, 4066, 342, 5, 24, "Output"],
Cell[156055, 4073, 563, 12, 18, "Message"],
Cell[156621, 4087, 561, 12, 18, "Message"],
Cell[157185, 4101, 561, 12, 18, "Message"],
Cell[157749, 4115, 578, 12, 18, "Message"]
}, Open  ]],
Cell[CellGroupData[{
Cell[158364, 4132, 141, 3, 24, "Input"],
Cell[158508, 4137, 74, 1, 24, "Output"]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)
